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Abstract: We consider the problem of soft gluon resummation for gauge theory ampH- 
tudes and cross sections, at next-to-eikonal order, using a Feynman diagram approach. At 
the amplitude level, we prove exponentiation for the set of factorizable contributions, and 
construct effective Feynman rules which can be used to compute next-to-eikonal emissions 
directly in the logarithm of the amplitude, finding agreement with earlier results obtained 
using path-integral methods. For cross sections, we also consider sub-eikonal corrections to 
the phase space for multiple soft-gluon emissions, which contribute to next-to-eikonal log- 
arithms. To clarify the discussion, we examine a class of log(l — x) terms in the Drell-Yan 
cross-section up to two loops. Our results are the first steps towards a systematic general- 
ization of threshold resummations to next-to-leading power in the threshold expansion. 
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1. Introduction 

The eikonal approximation has played an important role throughout the history of gauge 
theory scattering amplitude calculations. It embodies the universal behavior of such am- 
plitudes in the limit where massless gauge boson momenta become soft, and it encodes 
the semiclassical nature of soft, long-wavelength radiation. Early applications to potential 
scattering and Regge theory (sec for example [1, 2, 3]) predate the standard model of 
particle physics. In the context of modern applications to QED and QCD, the eikonal 
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approximation is typically phrased in terms of Wilson lines: when radiated gauge bosons 
are soft, one may neglect the recoil of energetic particles; these can then be replaced by 
straight Wilson lines along the classical trajectories of the emitters (often called eikonal 
lines for this reason); this expresses the fact that soft gauge fields can affect energetic 
particles only by dressing them with a gauge phase. The fact that the eikonal approxima- 
tion can be recast in terms of expectation values of (non-local) operators such as Wilson 
lines points to a deep connection between the perturbative and non-perturbative aspects 
of gauge theory scattering, a fact that was recently highlighted by remarkable results in 
maximally superymmetric A/" = 4 Yang-Mills theory (for a review, see [4]). Recent signif- 
icant developments in the understanding of the all-order structure of infrared divergences 
in non-abelian gauge theories [5, 6, 7, 8] were also derived with the help of the special 
properties of the eikonal approximation. 

A crucial feature of the eikonal approximation is the fact that scattering amplitudes 
exponentiate, which means that it is possible to establish a set of simplified rules to compute 
the logarithm of the amplitude. As a consequence, low-order perturbative calculations can 
be employed to gain access to all-order information, which is of great interest both for 
theory and for phenomenology. In QED, for example, a scattering amplitude Ai dressed 
with multiple soft photon emission may be written [9] in the schematic form 



where A^o is the amplitude without soft photon radiation, and the sum in the exponent 
is over connected Feynman diagrams for soft photon emission. Exponentiation was later 
shown to hold also for non-abelian gauge theory amplitudes [10, 11, 12], where however 
the structure of the exponent is more complicated, due to the non-commutativity of color 
matrices. The non-abelian analogue of eq. (1.1) for soft gluon emissions from two energetic 
(hard) partons connected by a color singlet vertex is 



where the sum in the exponent involves Feynman diagrams W which are two-particle 
irreducible with respect to the hard emitting partons, and are called webs in the literature. 
The factors cw are modified color factors, which differ from the conventional color factors 
cw associated with the same diagrams through the standard Feynman rules. Recently, 
eq. (1.2) was extended [13, 14] to more general multiparton amplitudes. In that case 
webs do not allow for a simple topological characterization, but it remains true that the 
logarithm of the amplitude can be computed in terms of a subset of Feynman diagrams 
with modified color factors. 

The exponential structure of the amplitudes is a necessary ingredient for the resum- 
mation of soft gluon effects in many different high-energy cross sections. In general, in 
order to perform a rcsummation for a given observable, one must also require that the 
phase space for multiple soft gauge boson emission should factorize. If that happens, 
potentially large logarithms of ratios of physical scales, which often jeopardize the ap- 
plicability of perturbation theory, can be resummed in exponential form. Soft gluon re- 
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summation is well understood, and has been widely investigated using a variety of tech- 
niques [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25], leading to a vast array of phenomenological 
applications to many scattering processes in perturbative QCD. In this connection, it is 
important to recall that some of the large logarithms that arise in high-energy cross sec- 
tions, for example when scattering takes place near a partonic threshold, originate from the 
emission of hard collinear gluons, and thus cannot be captured by the soft approximation. 
In what follows we will try to clearly distinguish these two effects. 

The aim of this paper is to clarify the structure of next-to-eikonal (NE) contributions 
to scattering amplitudes (and eventually cross sections) . For an amplitude involving n soft 
gauge bosons with momenta ki {i = 1, . . . ,n), we imagine rescaling all soft momenta by 
a common factor tj. Then the eikonal approximation to the amplitude keeps the leading 
power in r/ as 77 —7- 0, which is typically 77^" for an amplitude involving n soft gauge 
bosons. The NE approximation organizes terms proportional to the next-to-leading power 
in this expansion, 7y^("^i). Essentially, one takes fej — )• for all but one gauge boson, whose 
momentum is kept to one power beyond the eikonal approximation. This eikonal expansion 
is related, but not identical, to the expansion of partonic cross sections near kinematic 
thresholds, which is used to organize Sudakov logarithms. Near partonic thresholds, one 
can typically parametrize the distance from threshold in terms of a single dimensionless 
variable ^, such that ^ — ;> is the threshold limit. Examples include ^ = W^/Q^ ~ (1 — .tbj) 
in DIS, ^ = 1 — Q'^/s in Drell-Yan, and ^ = 1 — T for, say, the thrust distribution in 
electron-positron annihilation. In all these cases, the partonic differential cross section 
receives contributions of the form 



We call the expansion in powers of ^ in eq. (1.3) a "threshold" expansion. In this terminol- 
ogy, NE corrections to scattering amplitudes start contributing at the level of bnm terms, 
that is at next-to leading order in the threshold expansion (NLT). On the other hand, in 
the presence of final-state jets of hard particles, there are hard collinear emissions in the 
matrix elements which contribute at NLT but are not captured by the eikonal expansion: 
they must be dealt with using collinear evolution equations. NLT contributions can be 
phenomenologically important (see for example [26]) and are conceptually interesting since 
they probe the reach of the universal properties of soft and collinear radiation. As a conse- 
quence, these contributions have been recently studied in some detail: Ref. [27] showed how 
^-independent terms completely exponentiate for simple color-singlet QCD cross-sections 
(see also [28, 29]); subsequently, Ref. [30] proposed a modification of collinear evolution 
equations that allows the exponentiation of a subset of NLT terms; this was exploited in 
Ref. [31] to propose an ansatz for NLT Sudakov resummation for color-singlet QCD cross 
sections, that was shown to reproduce the bulk of these contributions through NNLO. 
Other studies [32, 33, 34, 35, 36, 37, 38] have uncovered intriguing relations between the 
coefficients of NLT contributions in finite order calculations, and have proposed all-order 
generalizations. The ultimate goal of our work, as well as of Refs. [39, 40], is to reach a 
complete understanding of NLT terms and organize them in resummed form to the ex- 
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tent to which this is possible. For cross sections, this wih also require an understanding 
of multi-gluon phase space beyond the eikonal approximation, and some progress in this 
direction is described already in this paper. 

As a first step in this program, we intend to establish to what extent the exponentia- 
tion properties of the eikonal approximation extend to NE level, and eventually to derive 
effective Feynman rules and iterative methods to compute directly the perturbative expo- 
nent to the required accuracy. Our results will reproduce and clarify from a diagrammatic 
viewpoint those obtained in [39], where the problem of NE contributions was first consid- 
ered using path integral methods. In that paper, a factorised form was assumed for Green 
functions with a fixed number of hard outgoing particles, which may emit any number of 
soft gluons. By recasting the propagators for the external particles (in the background of 
a soft gauge field) in terms of first-quantised path integrals, an effective field theory for 
the soft gauge field arises, with source vertices localised along the external lines. Within 
that formalism, exponentiation of soft photon contributions (in an abelian field theory) 
was shown to be straightforwardly related to the exponentiation of connected diagrams in 
quantum field theory. In the non-abelian case, the field theory for the soft gauge field was 
more complicated, due to the fact that the source terms arc matrix-valued in color space 
and, thus, non-commuting. It was, however, possible to ascertain that a subset of diagrams 
(indeed the webs of [12]) docs exponentiate. The exponentiation derived in [39] leaves out 
those contributions to the matrix elements that violate the assumed factorization prop- 
erties of the correlators. We similarly observe that our analysis applies to both real and 
virtual gluons, provided all emissions can be considered soft with respect to hard scale of 
the problem; in the case of hard virtual gluons, however, the known factorization properties 
of soft radiation do not fully generalize to NE level, and there are further contributions 
arising from soft gluons emitted by lines internal to the hard subamplitude. Here we will 
concentrate on the structure of factorizable contributions, which do exponentiate. 

To be specific, the diagrammatic techniques developed in the present paper, in agree- 
ment with the final result of [39] , organize next-to-eikonal contributions to matrix elements 
in the schematic form 

M = Mo exp [Meik + A^ne] (1 + Mr) + (NNE) . (1.4) 

Here A^o is the Born contribution, and Meik,NE collect factorizable contributions due to 
emissions of soft gluons external to the hard interaction. Non-factorizable contributions, 
arising from internal emission graphs, are collected in the function Mr- This contribution 
does not formally exponentiate, but has an iterative structure to all orders in perturbation 
theory. At NE level, as suggested in [39] , non-factorizable contributions can be organized by 
means of the Low-Burnett-Kroll theorem [41, 42], as extended in [43] to encompass collinear 
divergences. In the present paper, we neglect non-factorizable contrbutions contained in 
Mr and concentrate on the structure of the exponent in eq. (1.4) at NE level. Note that 
Ref. [39] addressed explicitly the case of two external hard partons connected by a color- 
singlet interaction. The structure of the exponent in the case of multi-parton amplitudes 
was recently examined in [13, 14], including a first analysis of next-to-eikonal corrections 
using path integral methods in [13]. 
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We emphasize that the diagrammatic analysis presented in this paper is a necessary 
step, not only in order to test the results derived with path-integral methods in [39], but 
also in order to place these results in the context of general proofs of factorization theo- 
rems, and to clarify the precise limits of applicability of the formalism. As an example, the 
factorization arguments organizing Low's theorem contributions at NE level in Ref. [43] 
are formulated in diagrammatic language, and a proper matching with factorizable con- 
tributions in order to avoid double counting must also employ diagrammatic techniques. 
Similarly, possible extensions of this formalism to collinear contributions at next-to-leading 
order in the threshold expansion will need to be mapped to existing diagrammatic anal- 
yses. Our results thus strengthen the validity of the conclusions reached by the path 
integral method, and place them in the more general context of factorization studies with 
diagrammatic methods. 

In order to test our results and to illustrate and clarify our technique in a concrete 
application, in the last part of this paper we consider a subset of the real-emission contribu- 
tions to the Drell-Yan cross section up to NNLO. Processes involving parton annihilation 
into an electroweak final state, such as Drell-Yan or weak boson production, or Higgs pro- 
duction in the gluon fusion channel, are ideal testing grounds for our formalism: indeed, 
near threshold all radiation in these processes is soft, so that one does not need to worry 
about hard collinear emissions. Furthermore, by considering real emission diagrams only, 
we don't need to include non-factorizable contributions related to Low's theorem, and 
therefore we should be able to reproduce exactly the fixed-order results upon applying our 
effective Feynman rules at the level of the exponent. We perform this test on abelian-like 
contributions (proportional to Cp in QCD) up to two loops (the first non-trivial order 
where exponentiation has an impact), finding the expected result. Applying our formalism 
to a hadronic cross section forces us to tackle the issue of the factorizability of soft gluon 
phase space, and we present some preliminary results suggesting that the factorization 
properties of the eikonal phase space may admit a simple generalization at NE level. 

The structure of the paper is as follows. In Sec. 2 we review the derivation of the 
exponentiation of soft corrections in the abelian and non-abelian cases, using an iterative 
Feynman diagram method. In Sec. 3 we derive effective Feynman rules for NE emissions. 
Crucial to this result is the fact that the sum over all NE diagrams has a factorisable form. 
In Sec. 4 we then demonstrate the exponentiation of NE corrections, extending the methods 
used in Sec. 2. In Sec. 5 we perform a detailed comparison of the rules thus obtained with 
the results found using the path integral method of [39]. Finally, in Sec. 6 we present 
the application of the ideas derived in this paper to the case of Drell-Yan production. In 
Sec. 7 we briefly summarize our results, while some technical aspects are discussed in the 
Appendices: notably. Appendix B discusses the structure of phase space for multiple soft 
gluon emission at arbitrary order, providing arguments for a form of factorization that 
could be applied at NE level. 

2. Introduction to eikonal exponentiation 

Exponentiation of soft gauge boson corrections at the eikonal level was first demonstrated 
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Figure 1: Soft photon emission from an energetic line. 



for abelian gauge theories in [9], and later generalised to non-abelian theories in [10, 11, 12]. 
Here we review the derivation of this result using Feynman diagram methods (see also [44] 
for a pedagogical exposition), both in order to make our paper reasonably self-contained 
and to introduce methods and notations that will prove useful when generalising the results 
to NE order. 

The proof of exponentiation is two-fold. First, one establishes a set of effective Feynman 
rules in the eikonal approximation, showing that they capture the leading-power behavior 
of the amplitude as soft momenta become vanishingly small, to all orders in perturbation 
theory; then one uses these effective Feynman rules to classify a subset of Feynman diagrams 
that generate the full eikonal amplitude upon exponentiation. We will do this first for the 
simple case of an abelian theory, and then discuss the non-trivial modifications that are 
necessary in the nonabelian case. 

2.1 Abelian eikonal exponentiation 

In order to derive the eikonal Feynman rules it is sufficient to consider a single hard 
massless external line of final on-shell momentum p, originating from some unspecified 

hard interaction described by a matrix element A^o(p)- The hard line may emit a number 
n of soft photons with momenta ki, as depicted in fig. 1, where we choose the ordering so 
that momentum ki is emitted closest to the hard interaction. If the emitting particle is a 
Lorentz scalar, such an external line dresses the hard interaction according to 

M^'-^-ip, h) = Moip) , A^2 (2f^ + ^1 + ^2)^^^ . . . r^V^(2p + Kr^r- , (2-1) 

[p + Kiy [p + Kny 

where we have introduced the partial momentum sums = Ylm=i ^'m- 

The eikonal approximation in this case can simply be defined as the leading-power 
contribution to the amplitude when the photon momenta A;^' — > 0, Vi. In this limit, 
eq. (2.1) becomes 

M^--(p,k)-MM ^^^-l^^^ . (2.2) 
If the emitter is a Dirac fermion, eq. (2.1) becomes 

M^-^-{p, h) = Mo{p) IIkIy ^'"' ■ ■ ■ ^pIkIy ^'"<p) ' (2.3) 

where u{p) is the spinor associated with the final state on-shell particle. At leading power 
in the soft momenta, this reduces to 

M>''-''-{p, h) = Moip) 7rA^7^^ • • • ir^l^'-uip). (2.4) 

2p ■ Ki 2p ■ J^n 



This appears to be more complicated than the scalar case due to the non-trivial spinor 
structure. One may however use the anticommutation properties of Dirac matrices and the 
Dirac equation to reduce eq. (2.4) precisely to the form of eq. (2.2), with the spinor u{p) 
reabsorbed into the radiationless matrix element M.o{p). The eikonal factor is the same in 
both cases, displaying the well-known result that the eikonal approximation is insensitive 
to the spin of the emitting particles. One may also notice that the eikonal factor docs 
not depend on the energy of the emitter, since it is invariant under rescalings of the hard 
momentum p^: at leading power in the soft momenta, one is effectively neglecting the recoil 
of the hard particle against soft radiation. 

The eikonal factor can be further simplified by employing Bose symmetry. Indeed, in 
constructing any physical quantity depending on the amplitude A^^^ "^" (p, A;;), one will 
need to sum over all diagrams corresponding to permutations of the emitted photons along 
the hard line. Having done this, the eikonal factor multiplying Moip) on the r.h.s. of 
eq. (2.2) will be replaced by the symmetrized expression 

£;«-^"(p,A;,) = -p^^...p^^- \, , ■■■ J ^, (2.5) 

TTl ~r fCjj-jJ P ■ {i^Tri "T • • • T~ i^TTn) 

where the sum is over all permutations of the photon momenta, and k^^^ is the i^^ momentum 
in a given permutation. There are n! permutations, and each gives the same contribution 
to any physical observable. This becomes manifest using the eikonal identity 

^1 1 1 -JJ ^ (2 6) 

Using eq. (2.6), the eikonal factor Ef^'^-^^"{p, ki) arising from n soft emissions on an external 
hard line becomes simply 

E''^-^-ip,ki) = (2-7) 

which is manifestly Bose symmetric and invariant under rescalings of the momenta {pi}- 
In practice, each eikonal emission can be expressed by the effective Feynman rule 



p-k 



As is well known, these Feynman rules can be obtained by replacing the hard external 
line with a Wilson line along the classical trajectory of the charged particle. In abelian 
quantum field theories this is given by the operator 



$^(0, oo) = exp 



Jo 



dX p ■ A{\^) 



(2.9) 



where /3 is the dimensionless four-velocity corresponding to the momentum p. This ex- 
presses the fact that soft emissions affect the hard particle only by dressing it with a gauge 
phase. 
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Figure 2: A process involving two eikonal lines A and B, interacting through the exchange of soft 
gluons forming diagram G. 



Having constructed the effective Feynman rules, one may proceed to demonstrate the 
exponentiation of soft photon corrections as follows. As an example, we consider graphs 
of the form shown in fig. 2, at a fixed order in the perturbative expansion. Fig. 2 consists 
of two eikonal lines, labelled A and each of which emits a number of soft photons. 
Notice that although we consider only two eikonal lines here, the analysis below generalises 
to any number of hard charged particles. One may envisage lines A and B as emerging 
from a hard interaction, and one may consider the graph G either as a contribution to an 
amplitude, or to a squared amplitude (in which case some of the propagators in G will be 
cut). 

Diagram G can be taken as consisting only of soft photons and fermion loops: in fact, at 
leading power in the soft momenta, soft photons originating from the hard scattering give no 
contribution (a statement which will have to be revisited when including NE corrections). 
Photons originating from one of the two eikonal lines must land on the other one, or on a 
fermion loop inside G. Indeed, a photon cannot land on the same eikonal line, as in that 
case the diagram is proportional to p^p^ = 0. 

Using eikonal Feynman rules, one finds that graphs of the form of fig. 2 contribute to 
the corresponding (squared) amplitude a factor 



where ki, Ij are the momenta of the photons emitted from lines A and B respectively, with 
i = 1, . . . ,n and j = 1, . . . ,m. 

Given that we have already summed over permutations in order to obtain the eikonal 
Feynman rules, each diagram G can be uniquely specified by the set of connected sub- 
diagrams it contains, as indicated schematically in fig. 3, where each possible connected 

(i) 

subdiagram Gi: occurs Ni times. According to the standard rules of perturbation theory, 
diagram G has a symmetry factor corresponding to the number of permutations of internal 
lines which leave the diagram invariant. This symmetry factor is given by 




G. 



>l.../in;t^l---t^m 




(2.10) 



(2.11) 
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Figure 3: Decomposition of a soft photon graph into connected subdiagrams Gc , each of which 
occurs Ni times. 



(i) 

where Si is the symmetry factor associated with each connected subdiagram or, and 
the factorials account for permutations of identical connected subdiagrams, which must be 
divided out. Contracting Lorentz indices as in eq. (2.10), the eikonal factor Tab may be 
written as 



= E n 1^ ki" 



NA 



Ni 



where 




(2.12) 



(2.13) 



is the expression for each connected subdiagram, including the appropriate symmetry fac- 
tor. Recognising eq. (2.12) as an exponential series, it follows that 



J^AB 



exp 



(2.14) 



We conclude that soft photon corrections exponentiate in the eikonal approximation, and 

the exponent is given by the sum of all connected subdiagrams. 

The same result was recently rederived in [39], using path integral methods to de- 
rive soft photon exponentiation from the well-understood exponentiation of disconnected 
diagrams in quantum field theory. 

2.2 Non-abelian eikonal exponentiation 

We now consider the generalisation of eikonal exponentiation to non-abelian theories. The 
results are well known [11, 12], but we will introduce methods and notation that will 
prove useful in the extension to NE order. Wc will examine the case of two incoming or 
outgoing hard emitting particles connected by a color-singlet hard interaction, as happens 
for example in Drell-Yan production, deep inelastic scattering or e+e" annihilation. We 
note that the results below extend to multiple hard colored emitters, as recently shown in 
Refs. [13, 14]. 

The proof of the abelian result relies crucially on the application of the eikonal identity 
after summing over the permutations of all photon momenta on each eikonal line. In the 
nonabelian case, this identity cannot be used, due to the presence of non-commutative 
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Figure 4: Decomposition of a higlier order non-Abelian diagram in terms of webs. 
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Figure 5: Gluon groups obtained from tlie lower eikonal line of fig. 4, with the same color coding 
as in that figure. 



color matrices associated with each emission. Exponentiation, however, is still possible, 
but with a somewhat more complicated structure. 

First we introduce the concepts of webs and groups. A web is a two-eikonal irreducible 
diagram, i.e. a diagram that cannot be disconnected by cutting the two eikonal lines^. 
Higher order diagrams can be rewritten as sums of products of webs: an example is illus- 
trated in fig. 4, which shows a particular fourth order diagram that is not a web, and its 
subsequent decomposition into webs. 

A group is the projection of a web onto a single eikonal line: therefore, gluons emitted 
from an eikonal line belong to the same group if they belong to the same web. This is 
illustrated in fig. 5, where we depict the groups that result on the lower eikonal line from 
the webs of fig. 4. 

Having introduced the concept of a group, one observes a useful generalisation of the 
eikonal identity, eq. (2.6), to sums over permutations that do not affect the ordering of 
gluons within groups [12, 44]. Let vr be a permutation of the gluons (with momenta {ki}) 
emitted by a given eikonal line of momentum p, with the restriction that the ordering of 
the gluons in each group g be held fixed. Then one may write 



I ... I = n 



p ■ ^Tri P • (^TTi + ^772) P ' {k^i + • • • + kff^) groups g P ' P ' (^91 ^92^' 



p ■ {kg^ + . . . + A;, 



(2.15) 



Here fe^-. and kg^ are the momenta of the z*^ gluon in permutation vf and group g, respec- 
tively, where g contains, say, m gluons. As an example of this result, consider an eikonal 
line with 3 gluon emissions, where gluons 1 and 2 belong to one group, and gluon 3 to 



^As shown in Refs. [13, 14], this simple topological identification of webs must be suitably modified when 
more than two eikonal lines are present. 
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another one; then, eq. (2.15) amounts to the statement 

11 1 11 1 

p- kip- {ki + k2)p- {ki + ^2 + ^3) p-kip- {ki + k^) p ■ {ki + /cs + ^2) 
_^1 1 1 1 

p- k2,p- {k^ + kijp- {ks + ki + k2) p-k^ \p-k2p- {ki + ^2) 

The RHS displays the factorization of gluons from different groups. Henceforth, unless 
otherwise stated, all permutations of gluon momenta involve fixed orderings of the gluons 
inside each group; thus we drop the tilde on the permutation symbols tt for brevity. 

Next we examine contractions of soft gluons emanating from two emitting lines con- 
nected by a color singlet hard interaction. A given diagram G connecting the external 
lines then has two color indices in some representation of the non-abelian gauge group, one 
index for each eikonal line. Color conservation forces the color structure of each diagram 
to be proportional to the identity matrix 5ij, where i and j are indices in the chosen rep- 
resentation. Denote by E(G) the eikonal amplitude arising from diagram G, not including 
color matrices. E(G) encodes the momentum information carried by a particular diagram, 
but not its color structure: we can then write the factor contributed by G to the (squared) 
matrix element as 

Tg = caSijEiG), (2.17) 

where cq is the color factor associated with G, and for simplicity we are not displaying the 
internal Lorentz structure. This notation allows us to rewrite eq. (2.15) in a more formal 
and useful way. Consider two soft gluon diagrams Hi and H2, connecting the same two 
eikonal lines A and B. One may define the product of these two eikonal diagrams, using 
eq. (2.15) in reverse, as 

E{Hi) E{H2) = E E H2) , (2.18) 

where Hi H2 denotes the particular merging that results from combining Hi and H2 
so that the gluon permutations on the eikonal lines are given by tta and ttb respectively^. 
Note that there may be a different number of gluons on each external line. An example of 
this is shown in fig. 6. 

Each possible merging of Hi and H2 gives some diagram G. The same diagram G 
could however result from different mergings that arise from different permutations tta, 
ttb- One may then rewrite eq. (2.18) more formally as 

E{Hi)E{H2) = E(G) Ngih.h, , (2-19) 

G 

where the multiplicity Nq^h^h^ denotes the number of ways in which diagram G can be 
generated from mergings of Hi and H2. Equation (2.19) generalises immediately to the 
product of any number of eikonal subdiagrams, as 

E{Hiy^ ■ ■ ■ mnY" = E ^ic) ^Vcii^r i/,^^...//- > (2-20) 

G 

^This product endows the set of eikonal gluon amplitudes with the structure of a shuffle algebra (for a 
definition, see for example Ref. [45]). 
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Figure 6: Illustration of the merging product UJ^ between two eikonal graphs, where tta and ttb 
are permutations of the gluons on eikonal lines A and B such that ordering within each group is 
held fixed. 

where each diagram Hi occurs Sj times in the product. 

With these notations in hand, we are now ready to state the exponentiation theorem 
for nonabelian eikonal diagrams. Consider the exponential 




(2.21) 



where the sum is over all diagrams H, each with an accompanying color factor ch whose 
interpretation will become clear in what follows (the bar is used to distinguish these color 
factors from those arising from conventional Feynman rules, {ch})- Note that each H may 
be decomposable in terms of products of smaller subdiagrams. Using eq. (2.20), one may 
then write 

exp|j;c-HE(i7)| =Ii(y.1[chE{H)A = J;cgE(G), (2.22) 

[ H J H \ n ^' / G 

where the sum on the right-hand side is again over all possible diagrams G {i.e. we use a 
different label on the left- and right-hand sides), and {cq} are constants. Expanding the 
exponential, one generates all possible products of diagrams. Each such product is itself 
equal to a linear combination of diagrams, by eq. (2.20), and thus contributes to various 
terms on the right-hand side of eq. (2.22). Each diagram H in the exponent thus appears 
on the right-hand side of eq. (2.22) in two ways: either as itself {i.e. contributing only 
to the term where G is equal to H), or as a component of larger diagrams G. One may 
then equate coefficients on the left- and right- hand sides of eq. (2.22), so that the choice 
of constants {cq} uniquely fixes, at least in principle, the constants {ch}- Our choice is 
to require that the constants {cg} on the right-hand side of eq. (2.22) be the usual color 
factors of perturbation theory. We can then show that the constants {ch} are zero except 
for a subset of diagrams H which have the property of being two-eikonal-line irreducible. 
These are the webs of [11, 12], that were referred to at the beginning of this section. 

One proceeds as follows. First one notes, as we just discussed, that each diagram G 
has a set of decompositions into subdiagrams H. This includes the case where G cannot 
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Figure 7: A diagram and its set of decompositions. The first decomposition has mi = 2 and 
m-x = 0, the second decomposition has mi = and mx = 1- 



be decomposed into smaller parts, and thus has only the trivial decomposition where H is 
equal to G itself. An example of how a simple diagram can be decomposed is shown in fig. 7. 
Each decomposition can be uniquely labelled by a set of integers {mn} that specify how 
many times each subdiagram H occurs in the decomposition of G. In the simple example 
shown in fig. 7, there are two decompositions labelled by {mi,mx} = {2,0} and {0,1}, 
where the notation mj^x is suggested by the figure. To make use of this decomposition, we 
note that one can rearrange the expansion of any product of exponentials, such as the one 
appearing in eq. (2.22), as 

n (e^ = E n^^ i'xm)V- (2.23) 

In the present context, each dummy index pn can be interpreted precisely as the multi- 
plicity ruH of subdiagram H in graph G. 

We can now exploit eq. (2.20) to rewrite each product of subgraphs in eq. (2.23) as a 
single sum over graphs. Equation (2.23) then becomes 

n fe [^^ E(i7)]«] = E(G) E ^o\in.,} f n --r] , (2-24) 

H \ n J G {run} \h ) 

where we denote by ^G\{mH} number of ways in which diagram G can be formed out of 
the given decomposition specified by {m//}. The expression now has the form of a single 
sum over all diagrams G (weighted by non-trivial coefficients), and thus can be matched 
to the right-hand side of eq. (2.22). One finds 




(2.25) 



This equation relates the coefficients {ch\ to the color factors {cg\- A similar relation was 
given in [44], and an explicit solution (in which the {c/fj are given in terms of the {cg}) 
was derived in [39]. One may interpret the coefficients \ch\ as modified color factors for 
the diagrams appearing in the exponent of eq. (2.22). At present this equation appears to 
contain no information, as the full set of subdiagrams appears in both the exponent and 
on the right-hand side. Crucially however, the modified color factors are zero except for a 
subclass of diagrams that are two-eikonal-line irreducible, as we now show. 

The proof proceeds inductively. The first step is to separate from the sum over de- 
compositions {m//} on the right-hand side of eq. (2.25) the term involving the trivial de- 
composition. Those remaining involve proper decompositions, where G genuinely reduces 
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Figure 8: General form of a diagram wliicli is two-eikonal-line reducible. 



to a product of smaller subdiagrams. This leads to writing eq. (2.25) as 

CG = -CG+Y: ^0\{m',} f n ^ > (2-26) 

where the prime denotes a proper decomposition, and the first term comes from the trivial 
decomposition. For example, for the case depicted in fig. 7, the proper decomposition is 
only the leftmost one {{m'j} = 2). Since cx = Cp — CfCa/'^; one finds cx = —CfCa/'^- 
One now assumes that the fact that the modified color factors are zero for two-eikonal-line 
reducible diagrams has already been shown up to some given order, so that lower-order 
diagrams appearing on the right-hand side can be taken to have this property. This will 
then be used to show that two-eikonal-line irreducibility persists to higher orders. 

Consider a general two-eikonal-line reducible diagram, shown in fig. 8. The normal 
color factor for such a diagram is given by cab = cacb- The sum over proper decomposi- 
tions for diagram AB has the form given in eq. (2.26), 



CAB = CACB = CAB + ^ ^ABlim 



where each subdiagram H is now two-eikonal-line irreducible by the induction hypothesis. 
This allows the multiplicity factor NAB\{m'jj} to be factorised. Each proper decomposition 
{m'jj} can be split into two parts denoted by {m^} and {m^}, where the subdiagrams 
in each case contribute solely to A or S respectively. That no subdiagram H contributes 
to both A and B follows from the two-eikonal-line irreducibility of the former. For each 
subdiagram H, one then has 

m'fj = m^ + mf[, (2.28) 

expressing the fact that the number of occurences of H in diagram AB is the sum of the 
number of occurrences of in A (m^) and the number of occurrences oi H in B (mf ). 
The number of ways in which H can occur in AB is then given by 

C(m^,mg)= a'^u B,r (2.29) 
(m^)!(mg)! 

which gives the number of ways of choosing occurrences of H, out of a total of m'jj. 
This allows one to rewrite eq. (2.27) as 



CA CB = CAB + -/V, 

{mi} {mf } 



Alimi} NBUm^^y (u ^ C(m^, mi) c™- j , (2.30) 
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where the original multiplicity factor -/VABI{m^} is replaced by separate multiplicity factors 
for A and B, times the number of ways of partitioning diagram H into A and B, for each 
H. Substituting eq. (2.29) into eq. (2.30), the double sum over proper decompositions 
appearing on the r.h.s of eq. (2.30) factorizes into a product of sums depending separately 
on A and B, and can be written as 




E 



B\{mB} 




(2.31) 



By eq. (2.25), this is equal to caCb, the product of the color factors of subdiagrams A and 
B. We conclude that cab = 0, which is the desired result. The fact that the lowest order 
diagrams are two-eikonal line irreducible completes the inductive proof. 

In this section we have derived the exponentiation of soft gluon corrections for two 
eikonal lines coupled by a color-singlet hard interaction, and in doing so have introduced 
methods and notations that will be useful in what follows. The nonabelian eikonal expo- 
nentiation theorem [11, 12] may then be summarised as follows. 

The sum of all diagrams involving soft-gluon emission from two eikonal lines 
connected in a color singlet exponentiates, and the exponent is a sum over two- 
eikonal line irreducible diagrams (webs, H) with modified color factors {ch}- 

In the rest of this paper, we consider the extension of this result to next-to-eikonal order, for 
the scattering of two hard partons. Arguments towards the generalization of these results 
to the case of several hard partons, using the results of Ref. [39], were given in Ref. [13]. 



3. NE corrections to Feynman diagrams 

In the previous section we reviewed the exponentiation of soft gluon contributions to scat- 
tering amplitudes. The proof of this result involves the use of effective Feynman rules in 
the eikonal approximation, valid at leading power in all gluon momentum components, as 
fc^' — 7- 0. The proof of the existence of effective Feynman rules in the eikonal approximation 
relies crucially on the factorization of contributions from individual photons (in the abelian 
case) and from different gluon groups (in the nonabelian case), expressed via the eikonal 
identity (2.6) and the generalised eikonal identity (2.15) respectively. 

We now wish to organize contributions to the amplitude that are subleading by one 
power of a soft momentum, which we call next-to-eikonal (NE) contributions. To be precise, 
as described in the introduction, if we rescale all gluon momentum components by a scaling 
parameter 77, as V^i, then the eikonal approximation to an amplitude with n (real or 

virtual) gluons retains only terms that are ©(r/^"), before loop momentum integrations. 
Here we wish to organize terms that are 0(77""+^). In order to do so, we seek to identify 
effective Feynman rules which can be applied in calculating diagrams where one extra 
power of gluon momentum is retained for at most one of the available gluons. We will 
see that it is indeed possible to obtain effective rules, and that this result again relies on 
the factorization of contributions from different gluon groups. The method used here is 
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Figure 9: The 2 gluon vertex arising in scalar field theories, which is absent for fermionic emitting 
particles. 

analagous to that used in the previous section, and will lead to a confirmation of the results 
obtained using the path integral approach in [39]. 

We will show that the effective Feynman rules for NE emissions include vertices for the 
correlated emission of gluon pairs, which correct for the incomplete factorization of indi- 
vidual gluon contributions at NE level. In a nonabelian theory, these vertices may couple 
gluons from different groups, and this will lead to the introduction of further contributions 
to the logarithm of the amplitude, arising from pairs of correlated webs. 

Beyond the eikonal approximation, it is no longer true that soft gauge boson emissions 
are insensitive to the spin of the emitting particle. For simplicity, we begin by studying 
the case of scalar particles. We will return to the spin 1/2 case in Sect. (3.2). 

3.1 NE emission from a scalar particle 

In the case of scalar emitting particles, the eikonal approximation to a scattering amplitude, 
given by eq. (2.2), receives NE corrections from three sources. 

• Subleading corrections to the gluon emission vertex, in which one of the momentum 
factors in the numerator of eq. (2.2) is replaced as 



• Taylor expansion of a propagator factor, which leads to the replacement of an eikonal 
propagator in eq. (2.2) as 



• For scalar emitting particles, the two-gluon ('seagull') vertex from the Lagrangian, 
shown in fig. 9, contributes to amplitudes at NE level. Indeed, a diagram with one 
seagull vertex has one propagator less then the corresponding diagram with only 
cubic vertices, so it is subleading by precisely one power of a soft momentum. 

The full NE amplitude consists of a sum over all possible insertions of the above replace- 
ments on each one of the hard lines. The seagull vertex has no analogue in the exact 
Feynman rules for emission from a fermion. We will see, however, that a vertex of this 
form arises in the effective Feynman rules for NE emissions from a fermion line, as shown 
in the following subsection. 






2{p ■ Ki)^ ■ 




p-Ki 
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In what follows, it will often be convenient to rewrite the NE emission vertex, using 
Ki + -fiTj+i = 2Ki — ki- Equation (3.1) then becomes 

p,, ^ Kl^i - ^kfi . (3.3) 

This decomposes the vertex into a part which depends only on a single gluon momentum, 
and a part that depends on a single partial sum of gluon momenta. Such a decomposition 
is useful, since any vertex which explicitly depends only on local gluon momenta can be 
immediately interpreted as an effective Feynman rule, as the other gluon emissions will 
factorize upon application of the eikonal identity, as in Sec. 2.2. 

To see this in more detail, note that a next-to-eikonal Feynman rule which depends on 
a single momentum in the numerator gives rise to a factor in the amplitude of the form 

^rn^. (3.4) 

j=l J 

where we have chosen the i^^ gluon to be NE, and we did not display the remaining 

numerator factors of p'^* [k ^ i) associated with the eikonal emissions. In evaluating the 
total NE contribution, by analogy with the eikonal case, one may replace the factor in 
eq. (3.4) with a sum over all permutations of the emitted gluons, holding the order in each 
group fixed. The total NE contribution the gives the factor 

^rEn^' (3-5) 

where, as in eq. (2.15), tt is a permutation in which the order of the gluons in a group is 
preserved. The product of denominator factors now factorises by the generalised eikonal 
identity, cq. (2.15), and one may extract the partial momentum sum Gi corresponding to 
the NE gluon (assuming this to come from group G) so that eq. (3.5) takes the form 

[...], (3.6) 

where the ellipsis denotes eikonal Feynman rules for all other gluons, factor ised into groups. 
The prefactor can then be interpreted as an effective Feynman rule. 

A similar argument applies to the two-gluon vertex of fig. 9. A diagram with such a 
Feynman rule contributes to the amplitude a factor 

n 

where again we are not displaying the factors of (A; / i, i + 1) from the eikonal emissions. 
One may sum over all gluon permutations tt, and the generalised eikonal identity then 
implies that the total contribution from the particular two-gluon vertex of eq. (3.7) gives 
the factor 

giJ-iiJ'i+i TT i (3 g\ 

groups 9 P'^aiP- (kgi + kg^) "' V (^91 + • ■ • + kg J ' 
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If the two gluons entering the two-vertex come from the same group G, one may extract 
the relevant denominator factors from eq. (3.8) to yield an expression of the form 



p-{Gi + Gi+i) 



where the ellipsis denotes a product of eikonal Feynman rules for all other gluons, factorised 
into groups. The prcfactor can be interpreted as an effective Feynman rule, coupling two 
gluons within a single group. If, on the other hand, the two gluons come from separate 
groups, then the seagull vertex couples together two groups: in other words, the denomina- 
tor of the effective Feynman rule contains the sum of two partial momentum sums Gi, Hj 
from different gluon groups. The two-gluon vertex therefore entangles two subdiagrams, 
each of which would be a web, were the seagull vertex not present. We return to this point 
in Sec. 4. Note also that the NE vertex still partly depends on sums over individual gluon 
momenta. This will be reconsidered in Sec. 3.3. 

3.2 NE emission from a spin-i particle 

In the case of spin-^ fermions, minimally coupled to massless gauge bosons, the lagrangian 
generates only three-point vertices. At NE order, we find therefore only two types of 
contributions, arising from corrections to the vertices and to the propagators respectively. 
Taylor expansion of the propagators in powers of the gluon momenta proceeds exactly as 
in the scalar case. The vertex corrections, on the other hand, are more complicated, owing 
to the non-trivial Lorentz and color structure. 

Let us first examine the spinor structure. Consider the numerator of eq. (2.3), where 
we drop the spinor u(p) and the hard matrix element Mo(p) for simplicity. We write 

M''''-'^-{p,ki) = {f+t^h^^^ ... {f+tn)!""^ =M^r''-{p)+N^]i-''-{pM) + --- , (3.9) 

where M^^'"'^" collects all terms linear in gluon momenta. In each such term, all factors 
of p to the right of the NE vertex insertion can be simplified using anticommutation and 
the Dirac equation (as in the eikonal case). We find then 



•'^NE 



{P,ki) = Yl 



j=i+i 



k=l 



(3.10) 



where we used anticommutation relations to move to the right of 7^^* in the last factor 
of eq. (3.10). The NE numerator is now written as the sum of two terms. The first term 
has no spinor structure, since one may anticommute all factors of through to the right, 
where they annihilate the Dirac spinor as in the eikonal case. One may then rewrite this 
term using 2Ki = Ki + i^i+i -|- fcj, as 



E 

i=l 



i=l 



lli2p^^)k^^ 



(3.11) 
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Figure 10: Schematic representation of the theorem given in eq. (3.12). Here • denotes a vertex 
which depends upon a sum of gluon momenta, and □ a vertex depending only upon the momenta 
entering the chosen vertex. 



The first summation on the right-hand side is the same as the sum over NE vertex correc- 
tions in the scalar case, as seen from eq. (3.1). The second summation is over terms which 
only depend upon a single gluon momentum, and therefore can be directly translated into 
Feynman rules. 

The second term in eq. (3.10), which still has a spinor structure, requires more careful 
handling. Using only the Dirac algebra, and a recursive argument, we prove in Appendix A 
that the next-to-eikonal numerator A/'^g'g '^"(p. A;,) can be written as 



i=l 



i-1 



E 

1=1 



(3.12) 



where the sum in the second term of the last line starts at i = 2. Equation (3.12) shows that 
spin-dependent terms in the sum over all possible NE insertions on a hard fermion line can 
be expressed as the sum of two kinds of contributions. The first kind implements single- 
gluon vertex corrections, each depending on the momentum of a single gluon. The second 
kind implements a new 2-gluon, 'seagull' vertex, as can be seen from the fact that the 
factor 2p ■ Ki in eq. (3.12) is precisely such as to cancel the eikonal propagator connecting 
the two neighbouring gluons i and i — The result is shown in diagrammatic form in 
fig. 10. 

Combining eqs. (3.11) and (3.12), one sees that the NE numerator for a fermion line 
is built out of two kinds of of contributions. There are terms where single and double NE 
emissions manifestly factorise {i.e. can be written in terms of the momenta and quantum 
numbers of the particles entering a specific vertex); there are however also terms which 
still depend on sums of gluon momenta along the eikonal line. We note that the latter, 
given by the first sum on the right-hand-side of eq. (3.11), together with the corrections 
from Taylor expansion of the propagators, are the same as in the scalar case, as seen from 
eqs. (3.1) and (3.2). 

Using arguments similar to those of Sec. 3.1, any term which depends only upon 
single momenta factorises according to the generalized eikonal identity, and thus can be 
interpreted as an effective next-to-eikonal Feynman rule. In particular, we find a one-gluon 
vertex resulting from the second sum in eq. (3.11) and from the first sum in eq. (3.12). 
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This vertex has the form 

V^^{k)^{k^-^^]^)T''. (3.13) 

In practice this will be combined with an eikonal denominator factor (involving the appro- 
priate partial momentum sum), as in eq. (3.6). Here, however, we focus on the numerator 
structure only, and postpone a full presentation of the effective Feynman rules to Sec. 4.1. 
We also find a two-gluon emission vertex, arising from the second sum in eq. (3.12), which 
has the form 

yj^^'«'' = 7M^^r"r\ (3.14) 

Note that we have explicitly reinstated the color factors associated with each emission 
in the case of a non-abelian theory. The contributions arising from eqs. (3.13) and (3.14) 
ultimately lead to a factoriscd expression for the amplitude. The same is not true, however, 
for the spin-independent contributions to the one-gluon emission vertex. We will see in the 
next section that the sum of these contributions can indeed be expressed in a factorised 
form, up to a remainder term expressing correlations between gluons in different groups. 

First, it is interesting to note some properties of the emission vertices in eqs. (3.13) 
and (3.14). Using 

l^^ = \{b^n + {l^M , (3-15) 
we see that eq. (3.13) may be rewritten as 

V^^ik) = -2ik,a^\ (3.16) 

where we recognise the generators of the Lorentz group a^^ = —\[y^,^'^]- The physical 
interpretation of this vertex, which does not occur for scalar emitting particles, is now 
clear: it is a (chromo-)magnetic moment vertex associated with the interaction between 
the spin of the emitting particle and the momentum of the emitted soft boson. This occurs 
for the first time at NE level, which is consistent with the fact that emitted radiation is 
insensitive to the spin of the emitter in the eikonal approximation. 

One may similarly replace eq. (3.14), using Bose symmetry, with its symmetric part 
under the exchange of the two gluons. One writes then 

The first term has the same form as the scalar seagull vertex, while the second term again 

involves the Lorentz generators a^^'^ . The latter term vanishes in abelian gauge theories, 
thus it contributes a non-abelian component to the chromo-magnetic moment. 

The results of this section can be translated into a set of effective Feynman rules de- 
scribing next-to-eikonal emissions. Some of these rules rely on individual momenta (they 
are 'local' in the pictorial representation of the Feynman graph), whereas others still rely 
upon partial sums of gluon momenta. Crucial to the eventual exponentiation of soft gauge 
boson corrections is the factorization of gluon or photon emissions. In the eikonal approx- 
imation, this factorization results from the application of the eikonal identity eq. (2.6) for 
abelian theories, and of the generalised eikonal identity eq. (2.15) for nonabelian theories. 
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We now consider the extension of these same methods to NE order. Having derived effec- 
tive Feynman rules, the next step is to show that contributions from different gluon groups 
factorise. As in the eikonal case, one may first separate the color structure of diagrams (as 
in eq. (2.17)), and then factorize the momentum structure arising from the NE Feynman 
rules. This factorization procedure is more complicated than at eikonal order (one lacks 
the relatively simple form of the generalised eikonal identity): in fact, it is necessary to 
introduce extra effective Feynman rules, which correlate the emissions of pairs of gluons 
from different groups. This is the subject of the following subsection. 

3.3 Factorization of next-to-eikonal diagrams 

In the previous section we have seen that there are two classes of effective Feynman rule 
at NE order: those depending only on the specific momentum associated with a particular 
gluon emission, and those that still rely on partial sums of gluon momenta. The first type of 
rule leads simply to the factorization of NE contributions from different gluon groups, using 
the same generalised eikonal identity, eq. (2.15), as in the eikonal approximation. More 
work is needed to analyze the contributions from the second type of Feynman rules. We will 
see in this section that factorization for these contributions is weakly broken by a remainder 
term, which implements correlations between pairs of gluon emissions from different groups. 
The resulting structure, however, is still sufficient to achieve exponentiation of NE effects 
in terms of a subset of diagrams, as shown in Sect. 4. 

Let us begin by introducing some notations and conventions that will be useful in the 
following. First of all, consistently with eq. (2.7), we denote the combined eikonal vertex 
and propagator by 

W{JP,K)^-^, (3.18) 

where one should keep in mind that in the following K will be allowed to be a partial sum 
of gluon momenta. 

Next, we introduce similar notations for NE corrections arising from numerators and 
from the Taylor expansion of propagator denominators. In accordance with eq. (3.1) and 
with eq. (3.2), we write 

NE^(p,X)^-^^-^p^ NE^(p,if)^-^, (3.19) 

where again K will generically denote a partial momentum sum rather than a single mo- 
mentum. One may also define the complete NE correction factor 

nW(jp, K) = NE^(p, K) + NE^(p, K) , (3.20) 

however, given that each diagram at NE order has either a vertex or a propagator correction 
(but not both), we will mostly consider the two types of correction separately in what 
follows. When considering subgraphs built with a single hard line, we will omit the hard 
momentum p as an argument in all these definitions. 

Consider next the organization of gluons into groups. As in the eikonal case, we 
partition the n gluons emitted by any one of the hard lines into groups g: gluons in the 
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Figure 11: Simple example of two groups of gluon emissions on a single eikonal line. See the 
discussion in the text. 

same group belong to the same two-eikonal irreducible subgraph (web) . Gluons are labelled 
along a given eikonal line as in fig. 1, i.e. from 1 to n, where the index increases as one 
moves away from the hard interaction. We will then need to sum over permutations vr such 
that the order of giuons in each group is held fixed. The order of gluons along the line in 
any such permutation is then given by vri, . . . ,7r,i, where tTj labels the i*'' gluon insertion 
moving away from the hard interaction, in permutation tt. 

It will also be useful to introduce a label I denoting those gluons which are the first of 
a group, in the order just stated. The group which contains this gluon we label by g{l). 
We then define 



In words, a tilde over a group denotes that group without the first gluon, if that gluon is 
I. Finally, as before, we find it useful to denote by G partial momentum sum of gluons 
in a group g (while we denote by K partial momentum sums in the entire set of available 
gluons). Thus Gk will be the partial momentum sum from gluon k up to gluon m, in a 
group g containing m gluons. We will also shorten Gi to G, representing the sum of all 
momenta in group g. 

To clarify the above definitions, consider the simple case of two gluon groups, depicted 
in fig. 11, where we label the various gluon emissions by upper and lower case letters, 
corresponding to each group. Fig. 11(a) shows a given hard line with two gluon groups 
g = {A, B, C} and h = {a, b}. In this permutation, one has tt = {A, a, B, C, h). Fig. 11(b) 

corresponds to a different permutation, vr = (^,S,a, 6, C), where the order of gluons in 
each group is not changed. In this simple case, the label / takes the values A or a, and 
g{A) = g, while g{a) = h. The partial momentum sums G and H in each case are given 
by G = kA + ks + kc and H = ka + fcfo. Finally, if I = A one has g{A) = {B,C} and 
h{A) = {a, b}; alternatively, if Z = a then one has g{a) = {A, B, C} and h{a) = {b}. 

In constructing a (squared) amplitude at NE level we will encounter expressions with 
the general structure of cq. (2.10), with two important differences: attachments of the 
gluon graphs to the hard lines will now include NE corrections, and, therefore, we will not 
be able to employ directly the generalized eikonal identity. As a consequence, it will be 




(3.21) 
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necessary to deal with sums over permutations, such as the one appearing in eq. (2.5). To 
handle these sums, we introduce a shorthand for the Lorentz tensors occurring in these 
expressions. For every allowed permutation it, we define 

n 

E(7r) = l[W-k{K^^), 



NE(7r) 



k=l 

n 

i=l 

n 



NE'Cvr) = ^ 



ki^i 



(3.22) 



NEi(7r) = NE'^-i (K^J n (K^J . 



k=2 

which will be convenient in the following proof. In this shorthand notation Lorentz indices 
are not displayed on the left-hand sides of eq. (3.22); also, one may, if needed, distinguish 
between propagator or vertex corrections for the NE terms, by simply adding the appro- 
priate index as in eq. (3.19). One easily sees that NE(7r) is the sum over all possible NE 
insertions, for a given permutation vr. Similarly, NE'(7r) is a sum over all insertions, where 
the first gluon is restricted to be eikonal, while NEi(7r) is the term resulting when the first 
gluon in permutation tt is next-to-eikonal. Thus one has 



NE(7r) = NEi(7r) + NE'(7r). 



(3.23) 



Finally, we note that we will employ the same notation as in eq. (3.22) for Lorentz tensors 
constructed with gluons belonging to a given group. Thiis, for example, E(^) will be the 
product of the eikonal factors associated with gluons belonging to group g (whose order is 
fixed in any permutation vr), while NE((7) will be the sum over all possible NE insertions 
in group g. Note that we also use the same notation for entire soft gluon diagrams, which 
were denoted by E{G) in the eikonal approximation already in Sec. 2, and will be denoted 
by NE(G) at NE level. This does not lead to ambiguities, since soft gluon diagrams are 
just products of line factors, such as the ones described here, contracted with suitable color 
tensors and other factors related to internal soft gluon interactions. These factors do not 
affect the manipulations performed here and below, which are related to the combinatorics 
of soft gluon insertions on the hard lines. 

We are now in a position to state the theorem governing the factorization of gluon 
groups at NE level. It can be written as 



Ene(^) = E 



NE(/i) n E(5) 



9+h 



R{g,h) W E(/) 



(3.24) 



This theorem can be seen as a further generalization, to NE level, of the generalized eikonal 
identity (2.15), which in the present language would read simply E(7r) = E(g'). The 
content of eq. (3.24) can be summarised as follows. 
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The NE contribution to soft gluon emissions from a hard line can be written in 
terms of two sums. In the first sum gluon groups are factorised, and each term 
contains a sum over NE gluon insertions in a given group, times a product 
of eikonal factors for the remaining groups. The second sum organizes pair- 
wise correlations between gluon groups, governed by a NE remainder function 
R{g, h), times a product of eikonal factors for the remaining groups. 

The proof is by induction, and it involves several steps. We must develop the induction 
argument separately for the NE vertex and propagator corrections: in each case, we will 
show a weak form of factorization which leaves behind pairwisc correlations of groups. The 
final step is to collect all these correlations and construct an explicit expression for the 
remainder function R{g, h), which will generate its own Feynman rule. 

We begin by noting that the theorem is trivially true for the emission of a single gluon. 
Next, we separate out the first gluon (the one emitted next to the hard interaction) from 
the sum over permutations. This gluon, being the first along the line, will also be the 
first gluon of a group. Thus, we can label it by I, and denote the group it belongs to by 
g{l). Furthermore, I will either be eikonal or next-to-eikonal. Then, using the notations 
introduced above, we may write 

^ NE(7r) = NE'^' (K) E E(^) + E (^) E ^^(^^ ' (^.25) 

TT I TT I n 

where K = ki is the sum of all emitted gluon momenta. In this equation, tt represents 
a permutation of all but gluon I, where the order in each group has been held fixed. Thus, 
we have split the sum over permutations into a sTim over the possible first gluons, and a 
sum over permutations of the remaining gluons. In the first term of cq. (3.25) one may use 
the generalised eikonal identity for the set of all gluons but the first, 

J2m=ll^{9)- (3.26) 

TT g 

In the second term of eq. (3.25) one may use the validity of the conjecture (3.24) for n — 1 
gluons. One then gets 

J2 NE(7r) = E NE« {K) J] m + ^ E'^' (K) I ^ NE(/.) J] E{g) + ...], (3.27) 

!■ 9 I \ h gj^h / 

where we did not explicitly write the remainder term, which will be dealt with in Sect. 3.4. 
In the next two subsections we consider separately the cases where NE''' represents the 
propagator and vertex corrections of eq. (3.19), and we prove that for both of these 
cq. (3.27) can indeed be reduced to the form of cq. (3.24). This will also implicitly determine 
the remainder term, which is then worked out in Sect. 3.4. 

3.3.1 The NE propagator 

In this subsection we consider the structure of eq. (3.27) with NE^(p,i^) NE'^{p,K), 
focusing on the propagator correction given in eq. (3.22). For brevity, below we will omit 
the dependence on the hard momentum p. 



- 24 - 



Let us begin by considering the first term in eq. (3.27). Here one may write 

2{p ■ K) 



^E^^{K) = ,^^-—^p.GiW^{Gi), (3.28) 



where the factor p ■ Gi cancels the eikonal propagator in E'*'(G/). Substituting this into 
the first term of eq. (3.27) and using 

E« (GO nE(5) = 11^(5)' (3.29) 

9 9 

one has 

J2 NE^' {K) n m) = (u Y.P-GI. (3.30) 

Now we may use the fact that Ylii Gi = K; indeed, note that Gi is the total momentum of 
gluons in group g{l), and the sum over I is equivalent to a sum over groups, since the order 
of emissions within a group is never rearranged. We can then write 

^ NE^' (K) n m) = J-J^ n • (3.31) 

I 9 ^9 



We may now use the fact that 



K^ = J2Gi + T.GrGr. (3.32) 

I lytr 



The cross-terms involve correlated emissions of gluons in different groups, and thus enter 

^2 



the remainder term, to be discussed in Sect. 3.4. For the Gf terms, one notes that 



NEi,p(5(/)) = ^^^^ P"' Em) = E{g{l)) , (3.33) 

where we absorbed an eikonal numerator and denominator in 'Ej{g{l)). One can use eq. (3.33) 
to obtain an expression for Gp which can then be substituted into eq. (3.31). One finds 

5;NE«(K) J{m = E 1^ NEi,p(5(0) n E(5) ■ (3.34) 

I g I ^ g^gil) 

Next, consider the second term in eq. (3.27). In this term, either h = g{l) or h ^ g{l), 
thus one may write 

J]E'^H^)ENEp(/i)nE(5) = EE'^'(^)NEp(5(0) n E(5) 

' h g^h l 9¥'9{l) 

+ NEp(/i) n E(5) , (3.35) 

I h^g{l) g^h 

where in the second term we used the fact that h does not contain gluon I to replace h 
with h. Next we use the simple identity 

E«(K) = ^E'^'(G;), (3.36) 

p ■ K 



-25- 



to find 

+ E E ^NEp(/.)nE(5). (3.37) 

In the first term on the right-hand side, we have used the definition in eq. (3.22) for 
NE'(g'(Z)) to recognise that 

NE'p(5(/)) = E« {Gi) NEp(^(Z)) , (3.38) 
while in the second term we have used the fact that 

E^^{Gi)llE{~g) = l[E{g). (3.39) 

Wc may now rccombinc cq. (3.37) with cq. (3.34). In doing so, the first term on the right- 
hand side of eq. (3.37) combines with the riglit-liand side of eq. (3.34) using eq. (3.23), and 
one finds 



5;NEp(7r) = ^^NEp(5(0) H ^(y) + E E NEp(/i) E(5) + • • • 

TT I ^ g^gil) I h^g{l) ^ g^h 



= E ^ E NEp(/i) n E(5) + • • ■ (3.40) 
I h g^h 

= ^NEp(/i) \{E{g) + ... , 

h g^h 

where the ellipsis again denotes the remainder term. In the second line, we have used the 
fact that = K . The result gives the contribution of the NE propagator corrections 

to the right-hand side of eq. (3.24), as expected. 

3.3.2 The NE vertex 

We turn now to the NE vertex correction NE^(p, K). Using 

K^ = G^ + Y,Gi, (3.41) 

we can rewrite the first term of eq. (3.27) as 

Y^^E^^{K)\[E{g) = Y,^-^^E,y{g{l)) \[ E{g) 
I g I ^ g^g{l) 

+ EEafnEe), (3.42) 

I k^l'- g 

where we have used the fact that 

NE^y{g{l)) = -±^E{~g{l))■ (3-43) 
P ■ 
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The second term in eq. (3.42) involves correlations between gluons in different groups, and 
enters the remainder term to be discussed in Sect. 3.4. The first term is of the same form 
as the analagous result for the propagator correction in eq. (3.34). 

For the second term in eq. (3.27), one proceeds as in Sect. 3.3.1. Eqs. (3.35 - 3.37) 
apply also in the case of the vertex correction, given that the form of NEp(5) is not used 
there. Given that the forms of eq. (3.34) and of the first term in eq. (3.42) are the same, 
combining the latter with cq. (3.37) leads (as in the previous subsection) to the right-hand 
side of eq. (3.24), which is the desired result. 

3.4 The remainder term 

We now turn to the determination of the explicit form of the two-group correlation function 
R{g,h) in eq. (3.24). As a first step, we note that we may again consider separately the 
propagator and vertex corrections, given that each NE diagram has at most one of these. 
Thus, one may write 

R{g,h) = Rpig,h) + Rv{g,h), (3.44) 

by analogy with eq. (3.20): the right-hand side contains the sum of the contributions due 

to propagator and vertex corrections respectively. 

Let us now rewrite eq. (3.27), including explicitly the contribution of the remainder 
term. One finds 

J2 NE(7r) = ^ NE« (K) J] ECg) + E^' {K) ^ NE(/i) ECg) 

I 9 I h gj^h 

+ ^ E'^' (K) R {g, h) n E(/) . (3.45) 

' 9^h fj^gA 

The derivation in Sects. 3.3.1 and 3.3.2 has shown that the terms in the first line 
of eq. (3.45) reconstruct the first sum on the right-hand side of our theorem, eq. (3.24), 
leaving behind two-group correlations which were dropped in eq. (3.32) and in eq. (3.41). 
Reinstating those terms, we see that what we have been able to prove so far takes the form 

^NE(^) = Y^mh) 11^(9) + Y^^'(^)Y.R{9,h) n E(/) 

^ h g^h I g^h f^g,h 

l¥=r ^9 l¥=r ^ g 

Separating the remainder function R{g, h) into its propagator and vertex components, 
according to eq. (3.44), we note that the first term in the last line of eq. (3.46) contributes 
(by definition) to Rp, whereas the second term contributes to Ry- Equating eqs. (3.46) 



-27- 



and (3.24) then gives the relations 

+ Ef^E(5)EW n E(5); (3.47) 

Y,Rv{g,h) n E{f) = ^E^^iK)Y,Rv{~g,h) U E(/) 

Qy^h f¥=g,h I g^h fy^g,h 

+ E^E(5)E(/.) n E(/). (3.48) 

which may be solved for the remainder terms Rpy. Notice that we have relabeled the sums 
in the second terms on the right-hand sides of eqs. (3.47) and (3.48) replacing 'Ylr^l ~^ 
^g^h, i-e. replacing the sum over all first giuons with a sum over all groups; also, we 
denoted by G and H the sum of all momenta in groups g and h respectively, and we 
introduced the notation g, which designates the group g without its first gluon (which 
implies that g = g if I £ g). 

We will proceed by writing down putative solutions to eqs. (3.47) and (3.48), and then 
showing that indeed the equations are satisfied. Our proposed solution to eq. (3.48) can 
be written as 

Let us explain the notations we have introduced. In the first sum, tt is any permutation of 
the giuons in groups g and h such that, as usual, the order of giuons in each group is held 

fixed. There is then a sum over all giuons in group g, where for each gluon one implements 
a NE emission vertex involving the partial momentum sum Hj for the gluon j from group 
h which lies nearest to the right of the gluon from g. That is 

where the individual gluon momenta hj all carry the Lorentz index of the gluon from g, 
so that the remainder term correlates this gluon with all the giuons in group h that lie 
to the right of the gluon from g. Finally, in eq. (3.49), E{ti) denotes a product of eikonal 
Feynman rules for all giuons except the gluon from g. 

To prove that eq. (3.49) solves eq. (3.48), let us begin by rewriting the latter as 

+ ^ E^' (K) ^ Rv (g, h) n E{f) , (3-51) 

' gy^h f^g,h 

In the second term, one may distinguish the cases in which gluon / is neither in g nor 
in h {I ^ {g, h}) from those in which I G g oi I G h. The second line in eq. (3.51) can 
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correspondingly be split into three separate terms, and one gets 

^RyM n E(/) -E^r^^^^f^EtejEW n E(/) 

g^h f^g,h g^h^ ^ ' i+g,h 



+ Y,Rv{g,h) ^ ^ n E(/) (3.52) 
g^h li{g,h} f¥^9,h 



+ ^^^-^^^^^^'^{G + H)Rv{g,h) n E(/) 

gj^h P f^g,h 

+ Y.^-^^-^^^''^G + H)^y^9,h) n E(/)- 

g^h ^ f^g,h 

Note that we have repeatedly used eq. (3.36) to rewrite E'^'- {K) in terms of group partial 
momentum sums. The fourth line contains all mergings of g and h such that the first gluon 
along the line (next to the hard interaction) comes from group h. This gluon does not 
couple to any of the gluons in g, since g and h arc assumed to be distinct; this is reflected 
in the fact that gluons from g couple only to gluons from h to their right-hand side, as 
described above. 

Similarly, the third line contains all mergings such that the first gluon along the line 
is a gluon from g, which does not coTiplc to any of the gluons in h. Note finally that, in 
the first line, the generalised eikonal identity implies 

ig,h) 

E(5)E(/i) = ^ E(7r), (3.53) 

TT 

where the right-hand side contains a sum over all mergings of g and h. The validity of 
eq. (3.49) can now be demonstrated (after some straightforward but tedious algebra) as 
follows. Substituting eq. (3.49) into eq. (3.52), the first and third lines combine to give all 
contributions where the first gluon along the line comes from group g. This then combines 
with the fourth line to give the sum over all possible mergings, with a prefactor in the sum 
oi p ■ {G + H)/p ■ K. Finally, this combines with the second term to give a total prefactor 

p-{G + H) ^ ^ p-G, 



p ■ K f-^ V - K 

Equating the result to the left-hand side of eq. (3.52), one finds that Rv{g, h) is given by 
eq. (3.49), which is then, as expected, the solution to eq. (3.48). 

We now proceed in the same way for the contribution of the propagator correction, 
Rp{g,h). Our proposed solution to eq. (3.47) is 



(9,h) 

^ 2p-{Gi + Hj)- 



Rp{g, h) = TyZ WZ^T^M^) ■ (3.55) 



Again there is a sum over all gluons in group g, where Gi is the partial momentum sum 
associated with this gluon. As in eq. (3.49), Hj is the partial momentum sum for the gluon 
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in group h which hes nearest to the right of the gluon from g. There is then a product of 
eikonal Feynman rules for all gluons in g and h, including those which are correlated by 
the two-gluon vertex, as is consistent with the first term in the last line of eq. (3.46). 
The proof is directly analagous to the vertex case. One first rewrites eq. (3.47) as 

g^h fj^g,h g^h > ^ f^g,h 

+ ^ E^' (K) ^ Rp{g, h) n E(/) , (3.56) 

and again considers the separate cases in which I G {g, K}, I G g and / G /i, which allows 
one to separate the second line of cq. (3.56) into three sums. Carrying out manipulations 
similar to those applied to eq. (3.52), substituting eq. (3.55) into the right-hand side, and 
combining terms, one finds as expected that Rp{g, h) on the left-hand side of eq. (3.56) is 
indeed given by eq. (3.55). 

This concludes our proof of eq. (3.24). We have shown that the sum over all possible 
NE gluon insertions exhibits a partial factorization into contributions arising from distinct 
groups. At NE level, a two-gluon vertex arises which correlates gluons at different positions 
along the line, including in general the case of non-adjacent gluons. This two-gluon vertex 
may either correlate pairs of gluons within the same group (in which case this remains a 
group), or it may correlate gluons in different groups (in which case the merging of g and 
h yields a single group). 

Having shown in this section that contributions from different groups factorise, one 
may proceed to show that next-to-eikonal corrections exponentiate. This is the subject of 
the following section. 



4. Exponentiation for NE matrix elements 

In Sec. 2.2 we showed that soft gluon corrections exponentiate in the eikonal approximation 
for non-abelian theories. Crucial to that derivation was the generalised eikonal identity, 
eq. (2.15), which states that, in summing over eikonal emissions, contributions from dif- 
ferent groups factorise on each external line. In this section we extend the argument to 
next-to-eikonal order, using the NE generalisation of eq. (2.15), given by eq. (3.24). 
By exponentiation at NE order, we mean the generalisation of eq. (2.22) to give 



^ [ciE(G) + cg^NE(G) 

G 



exp 



exp 



J^cj^(E(iJ)+NE(ii-)) 
L H 



H 



1 + ^ckNE{K) 



K 



(4.1) 



The left-hand side of eq. (4.1) consists of diagrams G, spanning the external lines, and 



evaluated up to NE order, accompanied by color factors c\ 



at eikonal (next-to- 



eikonal) order. On the right-hand side, one has an exponent involving diagrams H , where 
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the notations E(G) and NE(G) have the same meaning as in Sec. 2.2: they represent the 
momentum part of a given soft gluon diagram, with the hard interaction factored off. Care 
is needed in interpreting the sum over diagrams H. Those diagrams involving the next-to- 
eikonal single gluon and nonlocal two gluon vertices have topologies which are related to 
those already occuring at eikonal order (i.e. for the two-gluon case, we may think of this 
vertex as correlating two single gluon emissions). However, diagrams involving the local 
(seagull-like) two gluon vertex have topologies which have no counterpart at eikonal order. 
For such diagrams if, one has E(iJ) = in the above sum, so that the form of the above 
result is still correct. Note that eq. (4.1) implies that the modified colour factors for the NE 
diagrams are the same as those for the eikonal diagrams, when both give a non- vanishing 
contribution. We will see in what follows that this is indeed true. 

In the simple case we are considering (two hard partons coupled by a color singlet 
interaction), the color factors for all subdiagrams commute, since they must all be propor- 
tional to the identity matrix in the chosen color representation. As in the eikonal case, 
eq. (4.1) would contain no information if it were not for the fact that the modified color 
factors ch are zero except for a subset of diagrams, which are two-eikonal-line irreducible. 
The diagrams on the right-hand side can then be interpreted as eikonal and next-to-eikonal 
webs. Note that, up to NE order, one may rewrite the NE exponentiation theorem as shown 
in the second line, which consists of a Taylor expansion of the NE part of the exponential. 
Higher order terms in this expansion are NNE and so on, and thus one may choose whether 
to place NE webs in the exponent or not. We return to this point later. First, we must 
show that the form of eq. (4.1) is indeed correct. 

To simplify the argument, we begin by neglecting the remainder term in eq. (3.24). 
We also neglect the local two-gluon vertex discussed in sections 3.1 and 3.2, consisting of a 
seagull vertex, plus a spin-dependent correction for fermions. Thus, we consider only NE 
vertices involving the emission of a single gluon. By analogy with eq. (2.18), we may then 
use eq. (3.24) to write, for example, 

n¥.{HiH2) = NE(iii) E(ii2) + HHi) NE(ii2) = Y.Y^ ^^Hi H2) , (4.2) 

where NE(iJj) denotes the next-to-eikonal contribution to a subdiagram Hi (which is a 
sum over propagator and vertex corrections). We also introduce the notation NE(ifiif2) 
for the sum over all possible next-to-eikonal insertions in the diagram formed by a product 
of subdiagrams Hi and H2. As in Sec. 2.2, we use tt/ to represent a permutation of the 
gluons on eikonal line /, where the ordering of gluons in each group is held fixed. As stated 
above, we have ignored the remainder term R{Hi,H2) in eq. (4.2), so that only one-gluon 
vertices are present. Each merging of diagrams Hi and H2 leads to a new diagram G, but 
this diagram can be formed in different ways. As in the eikonal case, one may then write 

E{Hi) NE{H2) + NE{Hi) E{H2) = ^ NE(G) Ng\h^h, , (4.3) 

G 

where Nq^jj^jj^ is a multiplicity factor representing the number of ways that G can result 
from the merging of Hi and H2. Note that this is the same factor as in the eikonal 
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case, given that we are so far considering only one-gluon vertices, and the topology of the 
merged diagrams is independent of whether gluons are next-to-eikonal. Eq. (4.3) is easily 
generalised to the product of any number of diagrams Hi, where each diagram may occur 
ruHi times. Consider first the case in which one has mn copies of a single diagram H. 
Then one may write 

rriH-l 
rt=0 

= ruH NE(i?) E(i/)™«-^ , (4.4) 

where the notation NE(i^'"^) naturally generalizes the one introduced above for NE(i7ii72)- 
Here the first line follows directly from eq. (3.24) (in which the remainder term is also ne- 
glected), and in the second line we have used the fact that the factors E(//) and NE(i?) 
commute with each other. One may then write a general NE product of many diagrams, 
by analogy with eq. (2.20), as 



mu, NE(ffj) EiRj)""! JJ E(H,r'. 



j:mG)N^^^^„,^r.,^, (4.5) 

G 



where N u-^Hn is the multiplicity for forming diagram G from the set of sub- 

G\H-^ ...H„ 

diagrams Hi, each occurring mij. times. As discussed above, this is the same factor as in 
the eikonal case. 

Following the reasoning of Sec. 2.2, we now consider the expansion of the quantity 



exp 



H 



■ CK NE(i^) = Y,CK NE(K) • n I E [^hE{H)\ " I . (4.6) 

K K H K n ' J 



This corresponds to the second term on the right-hand side of eq. (4.1), and we will now 
show that the modified color factors ck indeed satisfy the same relation (in terms of the 
normal color factors ck) as in the eikonal case. This allows us to interpret the diagrams 
K as next-to-eikonal webs, which must also be two-eikonal line irreducible. 

In order to do this, we follow once again the reasoning of Sec. 2.2. Using eq. (2.23), 
we rewrite the product over diagrams H in eq. (4.6) as a sum over multiplicities mn- 
Then we examine a single term in the sum over K and {m//} that results on the r.h.s. of 
eq. (4.6). We rewrite this term by extracting from the product over diagrams H the factor 
corresponding to H = K. Wc get 

-cK^m)T{^,\cHm)r" = -cKmK)cr-' ^^^^""\] n ^^^^^^^ (4.7) 

mnl [piK - 1)! ^^j^ mni 

(\ -mj 
n E(F)™^ 
H^K J J 

where we also relabelled niK — > friK — 1, so that diagram K occurs a total of niK times. 
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Substituting this in eq. (4.6), and using eq. (4.5), one finds finally 



exp 



H 



^cxNE(K) = ^ NE(G) ^ iVG|{™} ( n 

G {ruH} 



K 



(4i 



H 



Equating terms of NE order on both sides of eq. (4.1), one finds that 



H 



(4.9) 



This relation between the normal color factors {cg} and the modified color factors {cg} 
is exactly the same as in the eikonal case, given by eq. (2.25). Thus, the modified color 
factors for NE webs are the same as those obtained for the corresponding eikonal webs. 
In particular, the modified color factors for two-eikonal-line reducible diagrams vanish also 
at NE level. This is not surprising, given that the color factor carried by a NE one-gluon 
vertex is the same as that carried by the eikonal vertex. The derivation of eq. (4.1), however, 
relies crucially on the non-trivial fact that the momentum-dependent contributions from 
different groups factor isc. 

Having demonstrated the argument for the simplest case in which correlations between 
groups (and the seagull vertex) are neglected, we now prove the result of eq. (4.1) in full. 
As shown in Sec. 3.4, the complete form of the NE generalised eikonal identity is given by 
eq. (3.24). The corresponding generalisation of eq. (4.5) is then 



+ E 

j 



mH,{mH,-l)R{Hj,Hj)mj) 



(4.10) 



mH.mH, R{Hi, Hj) B{Hir»i-^ E(iJ,) 



This formula includes three types of NE contribution: 



• Single gluon vertices: as before, these reside in the first term on the right-hand side, 
contributing to the factor NE(iJj). 

• The remainder term of Sec. 3.4: this contributes to the second and third terms on 
the right-hand side, via factors of R{Hj,Hi). 

• The local NE two- vertex of fig. 9 (in the scalar case) and cq. (3.12) (in the fermion 
case). As discussed in Sec. 3.1, this vertex may either couple gluons from the same 
group G, or it may couple gluons from two different groups G and H. In the former 
case, one includes it in the first term on the right-hand side of eq. (4.10); in the 

latter case, one regards it as part of the remainder term, since it couples together 
two groups, which would remain separate if the local two-vertex were not present. 
As remarked already above, E(if) = for such diagrams. 
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In those terms which involve the remainder function R, we expUcitly separated the con- 
tributions from correlations between two instances of the same subdiagram (involving 
R{Hj,Hj)), and correlations between different subdiagrams (involving R(Hi, Hj)). The 
factors in each sum come from the number of ways in which these combinations occur. 
Now, to proceed along the lines of the discussion leading to eq. (4.9), we consider the 
expansion of the quantity 



S = exp 



m,H 



H 



Ci^NE(iC) + RiK, L)cKCL \=£i+S2. (4.11) 



K 



K,L 



The first term in square brackets leads to eq. (4.7). The second term contains correlations 
between groups, and one may again separate out the term with K = L in the sum, to 
obtain 



£2 = 



K 



H K n 



(4.12) 



As above, one may rewrite the expanded exponential in terms of a sum over the possible 
multiplicities {mu}. In doing so, one may take the terms involving K and L out of the 
product over diagrams H. It is also convenient to relabel rriK and m^, where appropriate, 
so that each group H occurs mn times in each sum. Eq. (4.12) then becomes 



{rriH} 



+ 1^ckclR{K,L) ^ YV 



{rriK - l)\{mL - 1)! ^^^^^ mj/! 



. (4.13) 



Combining eq. (4.13) with eq. (4.7) gives 



{mH} 

+ Y,mK{mK -l)R{K,K)E{Kr^-^ n E{H) 
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K H+K 
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+ ^m;^mLE(i^,L)E(i^)"^^-iE(L)"^i-i J] E(H\ 



rriH 



K,L 



n 



mjl 



(4.14) 



The prefactor in brackets has precisely the form required by eq. (4.10), and one finds 



£ = exp 



rriH 



H 



Y CK NE(iC) + Y CKCL R{K, L) 



K 



K,L 



5^NE(G) Y ^G|{™} n 



{mH} 



H 



(4.15) 
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Comparing again with eq. (4.1), one finds precisely the same color relation already ex- 
pressed in eq. (4.9). This completes the proof of the relation 



as given in the second line on eq. (4.1). Here however, in a slight abuse of notation, we 
have explicitly extracted the remainder term in the ncxt-to-cikonal contribution on the 
right-hand side. The modified color factors in this equation are the same as those obtained 
in the eikonal case (for those topologies that occur already at eikonal order). As stated 
above, the NE terms can be moved into the exponent up to NE order. The content of this 
theorem can be summarised as follows. 

Matrix elements containing gluon emissions computed up to next-to- eikonal 
level exponentiate, and the exponent contains a sum over webs. The NE webs 
belong to two classes: (i) two-eikonal-line irreducible subgraphs including at 
most one NE emission vertex, and (ii) novel two-eikonal-line irreducible graphs 
resulting from joining together two eikonal webs with a two-gluon vertex. 

We emphasize again that the exponentiation theorem of eq. (4.1) applies only to contribu- 
tions from soft gluon emissions which are factorizable from the hard interaction. In order 
to gain complete control on NE logarithms in the presence of hard virtual gluons, one must 
still include the non-factorizablc corrections denoted by Mr in eq. (1.4). Furthermore, in 
order to gain control on the threshold expansion at NLO, one must also take into account 
hard coUinear gluon emissions, which lead to threshold logarithms of the same order as NE 
terms. We leave the study of these further contributions to future work. 

4.1 Next-to-eikonal Feynman rules 

In the preceding sections we have classified all next-to-cikonal contributions in diagrams 
for scattering amplitudes involving two colored lines, and shown that they exponentiate in 
terms of (next-to-)eikonal webs. In this section we provide a set of effective Feynman rules 
with which these webs can be efficiently computed. 

Our rules are defined elaborating upon those of the eikonal approximation. As usual, 
one begins by identifying eikonal webs through the criterion of two-eikonal-line irreducibil- 
ity. The momentum space part of each web is computed using the Feynman rules for 
the eikonal approximation, given in eq. (2.8), and the corresponding color factor can be 
computed either recursively [11, 12, 44] or directly, using the replica trick [39]. 

At next-to-eikonal order, three new types of webs appear. The first type consists of 
diagrams obtained by replacing precisely one vertex or propagator in an eikonal web with a 
next-to-eikonal one, which must be done for each vertex or propagator in the eikonal web. 
The second type is found by constructing new two-eikonal-line irreducible diagrams, using 
eikonal Feynman rules, and precisely one NE two-gluon vertex. Together, these two types 
of webs constitute the second term in the last factor in eq. (4.16). The third type of web 




(4.16) 
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corrsponds to the third term in the last factor in eq. (4.16). At order k, one must consider 
all pairs of eikonal webs of lower orders i,j such that i + j = k, and for each such pair 
implement the correlation in eqs. (3.49) and (3.55). The rules for the color factors of each 
next-to-eikonal type are given in the previous section. Let us now treat each type in more 
detail, giving the effective Feynman rules involved. Note that we will use a photon-like 
symbol to represent a generic gauge boson, which may be abelian or non-abelian. 

For the first type of webs, we can formulate a one-gluon vertex rule, combining the 
second term on the right-hand side of eq. (3.3) with the propagator correction of eq. (3.2). 
For scalar partons, the result is 



(4.17) 



Note that because we provide rules that are only to be used in web diagrams, NE insertions 
and the partial momentum sums that they contain are automatically restricted to a single 
group, unless otherwise stated. The spinor equivalent of eq. (4.17) is 

k 

- ^ \ 2fG, 2{p-G,rJ - 2fG,y^i ^'^i+l^ P-Gi ) ^ 



p+C, p+C,^ 



(4.18) 

which follows from eqs. (3.11) and (3.12). On the right-hand-side, we have shown explicitly 
how the vertex part of this Feynman rule can be thought of as the sum of a spin-zero piece 
and a magnetic moment vertex, as remarked in Sec. 3.2. 

For the second type of web we merely need to give the two-gluon vertex rules. Any 
two-eikonal-line irreducible diagram one builds with one such vertex and any number of 
eikonal vertices and propagators is a next-to-eikonal web. For the case of a scalar emitter, 
one has the NE vertex rule 



- 2p.(fc+0 ' 

P (4.19) 

where we have included the relevant eikonal propagator. Note that and are color 
generators in the appropriate representation and that the numerator corresponds in fact to 
the exact vertex of scalar perturbation theory, as given in fig. 9. For the equivalent spinor 
vertex one has 




- pik+l) ' 2pik+l) ' 

(4.20) 

where the second term is the chromomagnetic moment vertex referred to in Sec. 3.2. The 
momenta in eqs. (4.19) and (4.20) are single gluon momenta in the abelian case. In the 
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non-abelian case, one must use the appropriate partial momentum sums corresponding to 
each gluon: if the gluons both originate from group G (with numbers i and i + 1), then one 
should use k = Gi and I = Gj+i. If the gluons come from different groups G (gluon number 
i) and H (gluon number j), one should evaluate the Feynman rule using the appropriate 
partial momentum sums Gi and Hj. Note that in this case it is more natural to absorb 
this Feynman rule as part of the remainder vertex. 

Finally, in the third type of web we must correlate any two gluons in different groups, 
by means of an additional two-gluon vertex rule that corresponds to the remainder terms 
derived in Sec. 3.4. This has the form 



Here the tensor product in the color factor indicates the fact that the color matrices corre- 
spond to gluons on different parts of the external line, belonging to sets that in the eikonal 
approximation would be different groups. Recall that Gi and Hj denote the partial mo- 
mentum sums for gluons at position i (j) in group g (h), where it is understood that j > i, 
as described in Sec. 3.4. 

The above results match those obtained in [39]. It is useful to relate our results here 
with those of that paper where they were derived using path integral methods. This 
comparison is the subject of the following section. As it is not immediately relevant to the 
demonstration of the application of the Feynman rules in Sec. 6, the reader may omit the 
following section without interrupting the flow of the paper. 

5. Relation to the path integral formalism 

In this section, we briefly review the method and results of [39], which flrst demonstrated 
the exponentiation of a class of NE corrections using a path integral approach, and derived 
effective Feynman rules for NE emissions. A comparison of the two approaches is useful for 
two reasons. First, it offers a non-trivial check of the effective Feynman rules derived here 
and in [39]; second, the present derivation reveals certain subtleties associated with NE 
Feynman rules, that were not fully addressed in the path integral approach. For the sake 
of completeness, we will briefly describe the path integral approach, referring the reader 
to Ref. [39] for a full exposition. We shall moreover restrict ourselves to the case of scalar 
emitting particles, which will be sufficient to clarify the correspondence between the two 
approaches. 

The starting point in [39] was to separate the gauge field into hard and soft modes, 
and then consider (as wc do here) a scattering amplitude for the production of n final state 
hard particles, each of which may emit further soft radiation. For the simple case of abelian 
gauge theory with scalar emitting particles, such a scattering amplitude has the factorised 




(4.21) 
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form 

/„ n n 
P^^ / l[dxi H{xi, ...,xn) l[{pk\{S^u^ - isr'\xk)pl e^^[^^] . 
i=l k=l 

(5.1) 

Here A'^ is the soft gauge field (with action ^'[tIs]), and H{xi, . . . ,Xn) is the factor repre- 
senting the hard interaction producing the emitting scalar particles at positions Xi. The 
factors {pk\{Squa.d ~ ^^)~^\xk) represent propagators for a scalar particle in a soft back- 
ground gauge field (with Squad denoting the quadratic operator for the scalar field in the 
Lagrangian). These propagators are computed between states of given initial position (the 
point Xk at which the particle is created by the hard interaction) and given final momentum 
Pk- The explicit factors of p^. in eq. (5.1) truncate the free propagators associated with the 
external legs. As noted above, the factorization represented by eq. (5.1), which is exact in 
the eikonal approximation, does not encompass all NE corrections to the amplitude, since 
there are contributions arising from virtual diagrams containing hard gluons which do not 
factorize in the form of eq. (5.1). Such contributions can be traced and organized using an 
appropriate extension of Low's theorem [43] , and their treatment is left to future work. 

The propagator factors in eq. (5.1) can be represented as first-quantised path inte- 
grals [46, 47]. One finds 

pI (PikKVad - ^e)-'\xk) = e-'P'^-^^fkioo) , (5.2) 

where 

/jt(oo) = / Vyk expii / dt 

Jj/fc(0)=0 I Jo 

(5.3) 



^vlit) + {Pk + Vkit)) ■ As (xk+Pkt + yk{t)) 



+ ■ As{xk+ Pkt + yk{t)) 

Here yk{t) is the fluctuation about the classical straight-line path associated with the A;-th 
scalar particle, and the path integral is over all such fluctuations, subject to the boundary 
conditions of given initial position Xk and final momentum pk- The implementation of 
these boundary conditions is discussed in more detail in Ref. [39]. Substituting eq. (5.3) 
into eq. (5.1) yields 

Mipi,...,pn) = [VA>^ [ fldx, I f[PyfeiJ(xi,...,ar„)e^^[^^le-(-i-fi+-+-"-^'") 

•' i=l k=l 

X cxp i^ dt Qy^(0 + (pk + Vkit)) -As + ^d- As^ . (5.4) 

We see that the scattering amplitude has taken the form of a generating functional for a 
quantum field theory of the soft gauge field A^. The exponent in the second line contains 
terms linear in As, which act as sources. These sources are located on the external lines, 
and the path integral over Ag generates all diagrams that connect these sources. Thus 
formulated, the exponentiation of soft photon corrections (in terms of connected subdia- 
grams) is precisely equivalent to the well-known exponentiation of connected diagrams in 
quantum field theory [39]. 
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One advantage of this approach is the clear physical interpretation in terms of the 
worldline trajectories of the emitting hard particles. To see which diagrams exponentiate, 
one must calculate the soft gauge field Feynman rules that result after carrying out the 
path integrations over the fluctuations yk{t) in eq. (5.4). This was achieved in Ref. [39] 
by systematically expanding about the classical straight-line trajectory Xk{t) = + Pkt, 
for each external particle moving in the direction of p^. Note that the classical trajectory 
corresponds to the eikonal approximation in which the emitting particles do not recoil. 
Following Ref. [39], we rescale the hard momenta by means of a scaling variable A, so 
that, for large A, terms of order A° and A~^ constitute the eikonal and the next-to-eikonal 
approximations, respectively. With this choice, the expansion of eq. (5.4) up to NE order 
reads 



Mipi,...,pn) = jvA^J ndx,i/(xi,...,x„)e^^[^^le-(-i-^'i+-+-"-f") \[Fk{nk^As 



(5.5) 



where is the direction of the hard line carrying momentum p^, and 



dt 



Hk ■ As{nkt) + ■ As{nkt) + ^tnk^,UA'^{nkt) 



Fkink,As) = exp |i J 

/•oo /-oo r J ■ 
+ dtj^ dt' -—S{t-t')g^,A'^{nkt)A'',{nkt')-j9{t-t')nk^,g,Ad''A'^{nkt)] 



A'siukt')] - ^minMnk^nkug"^ [d^A'^inkt)] [a,A^(nfei')] 



(5.6) 



Eq. (5.5) now has the form of a path integral over As only, and the Feynman rules that 
couple the soft gauge field to the external lines can in principle be read off from the 
exponent. To that end, it is conventional to first perform a Fourier transform to momentum 
space. The practical procedure for this, in the present context, depends upon whether the 
gauge field is abelian or not. 

Consider first briefly the case where the gauge field is abelian, so that the source terms 
for As in eq. (5.5) commute. One may perform the Fourier transform in the exponent in 
the second line, as was done in Ref. [39]^. For a line with direction n^, one obtains 




+ 



(27r)"' J (27r) 
(fc • l)nf,nu 



d^l ~ ~ 1 
-A^{k)A'^,il)- 



n„ 
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+ V 



n ■ k X \2n ■ k 



n, 



2(n • ky 



rii/l^n • k + n^kyU • I 



\\2n-{k + l) 2{n ■ l){n ■ k)[n ■ {k + I)] 



2{n ■ l){n ■ k)[n ■ (fe + /)]. 



(5.7) 



with As{k) the soft gauge field in momentum space. The first line in eq. (5.7) contains one- 
photon emission vertices at both eikonal and next-to-eikonal accuracy. When combined. 



^Wc note some misprints in [39], involving the wrong correlators being used in eqs. (B.26) and (B.28) 
of that paper. The results as given here are correct. 
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Figure 12: A scries of emissions from an external line in position space. In a non-abclian theory, 
emissions are ordered, and path ordering of the color matrices is equivalent to time ordering of the 
emissions, with ti < t2 < ■ ■ ■ < tm- 



these correspond exactly to the rule of eq. (4.17). The second and third hnes of eq. (5.7) 
contain contributions to a two-photon vertex. Combining these terms results precisely in 
eq. (4.21), after removing color matrices and with the replacement Gi — t- k, Hj — t- I. 

Having seen how the correspondence between the two approaches works in the case of 
abelian gauge fields, we now consider the non-abelian case. As in Ref. [39], we assume a 
color singlet hard interaction with two outgoing lines. Performing the path integral over 
fluctuations yk{t), we get once again eq. (5.5), with n = 2. Now however the soft photon 
fields As in eq. (5.6) must be replaced with matrix valued soft gluon fields = As^aT"", 
which do not commute; as a consequence, the exponential in eq. (5.6) must be replaced 
with a path ordered exponential. 

To illustrate the consequences of this fact, consider for example fig. 12, which shows 
a number of eikonal emissions from the same external line, which would be generated by 
the first term in the exponent of eq. (5.6). For an abelian gauge field emitted from a hard 
particle moving in the direction n^, a line dressed as in fig. 12 would lead, after a Fourier 
transform to momentum space, to a factor in the amplitude of the form 



Leik{Tl, ki) — — r 

ml 



i J dim'^ie^^'^i*! ...i / dtmn^-^e 



OO 



{ki) ■ ■ ■ An^ {km) , 

(5.8) 

where henceforth for brevity we omit the subscript s on the soft gauge fields. The expo- 
nential factors in eq. (5.8) arise from the Fourier transform. Carrying out the U integrals'^, 
the coefficient of the product of photon fields (ki) . . . A^„(A;„) is 



1=1 ^ ' 



(5.9) 



which is the familiar product of eikonal Feynman rules for m photon emissions, with un- 
correlated momenta. Recall that in the diagrammatic approach the same result arises 
from the application of the eikonal identity, eq. (2.6), after relabelling photon momenta 
and summing over permutations. In the non-abelian case, the path ordering prescription 



''We implicitly assume the prescription n ■ ki + is to ensure convergence of the U integral. 
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modifies eq. (5.8) to yield 



Leik(n, fcj 



/"OO poo 

I / (iiin^ie''^-'=i*M / dt 
Jo Jti 

X A-m {ki) ■■■ A-nm (km) , 



(5.10) 



where the hmits of integration have been modified to reflect the time-ordering of the gluon 
fields. Carrying out the ti integrals, we find that the coefficient of A^^(A;i) . . . A^^(A;„) 
becomes 



(-1)"' 



{ki + ... + km) 



n 



n ■ Ki 



. (5.11) 



One observes that the effect of path-ordering is to replace single gluon momenta in the 
denominators of the eikonal Feynman rules with partial momentum sums. The sums are 
restricted to gluon groups since, as shown in Sec. 2.2, contributions from different groups of 
gluon emissions factorise. This indeed corresponds to the effective eikonal Feynman rules 
summarized in the previous section. 

Let us now analyse the remaining terms in the exponent of the non-abelian version 
of eq. (5.6) along similar lines, in order to recover the NE Feynman rules. Recall that 
the rules are different depending on whether gluons are in the same group or in different 
groups. Let us therefore first consider the case of gluon emissions along an external line (as 
in fig. 12), where the gluons all belong to the same group. We will see that different subsets 
of terms in the exponent of eq. (5.6) combine to give the Feynman rules obtained in the 
diagrammatic analysis of the previous sections. Let's begin by considering, in coordinate 
space, the second term on the first line and the second term on the second line of eq. (5.6), 
in the nonabelian case. They read 

/"OO r 1 

^'^E in,Xi) = / dt --d-A{nt) 

JO I 

poo poo 

-i / dt dt'e{t-t')ni'ga^[d''Ai^{nt)][A''{nt')] , (5.12) 
Jo Jo 

where we have set A = 1 in eq. (5.6). When the exponent is expanded (including the 
eikonal terms we discussed above), the first term in eq. (5.12) generates diagrams of the 
form shown in fig. 12, where one of the emitted gluons has one NE one-gluon vertex. If 
this is the j*'' gluon, the corresponding momentum space expression is 



-'NE 



{n, ki) 



poo / ■ \ poo poo 

I / diin'*ie'"-'=i*i ... — / dtjk''' e'^'-^i^i ...i 

Jo V 2/ Jtm-l 



dtmu'^'^e 



i.Tl'ituit'n 



X (fci ) . . . A^^ (fcj, 



(5.13) 



Carrying out the ti integrals, the coefficient of the product of gauge fields in eq. (5.13) 
becomes 

^NE,(1,J) ^^1,1^1) - ^ r. 11 



n- Kk 



(5.14) 
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Figure 13: Schematic form of diagrams involving a two-gluon vertex. Additional eikonal emissions 
may occur between the two correlated gluons. 



The second term in eq. (5.12), on the other hand, leads to diagrams of the type shown in 
fig. 13. In such diagrams, a pair of gluons is correlated through the appropriate two-gluon 
vertex. Note that the two gluons need not be adjacent on the external line, as indeed 

eq. (5.12) allows: path ordering restricts the second field merely to lie to the right of the 
first in the ordered gauge field product, and further eikonal emissions may lie in between. 
More specifically, taking the left-most correlated gluon to be the j'*^ one, fig. 13 leads to 
the momentum space expression 



i>j 



f'oo f'OO poo 

i / diin''ie''^-'=i*i ...I dtje'^'-^i^i . . . i / diin/^'fef e''^'^^'*' . . . 

Jo J tj—i J ti—1 



(A;i) . . . A^^ (km) ■ 



(5.15) 



After carrying out the ti integrals, the coefficient of the gluon fields in eq. (5.15) is 



n 



^>J l^j ^ ' 



n 



(5.16) 



The total contribution from the gluon, arising from the terms in eq. (5.12), is then 
given by the sum of eq. (5.14) and eq. (5.16), which gives 



^NE,(l+2,i) 



(n, ki) 



2n 



n 



n ■ Kk 



2n ■ Ki 



n 



n 



n ■ Kk 



(5.17) 



As expected, this result has the form of a product of eikonal factors times a NE Feynman 
rule for gluon j, and is precisely equal to the first part of the expression for effective NE 
Feynman rule in cq. (4.17). 

Next, we consider the last term in the first line and the term on the third line in the 
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exponent of eq. (5.6). In the non-abelian case, and for A = 1, they read 



E 



(3+4) 
NE 



1 f°° 

{n,Xi) = — / dttn''DAij,{nt) (5.18) 

2 Jo 

-- / dt dt'min{t,t')n>'n'' g"^ [daA^{nt)][drA^,{nt')] . 
2 Jo Jo 

The first term leads again to diagrams of the form of fig. 12, but with one of the eikonal 
emissions replaced by a NE one-gluon vertex, this time of the type generated by the ex- 
pansion of a propagator. Without loss of generality we can take this gluon to be thej*'* 
gluon. The corresponding momentum space expression is 



T (3j) 



(n, fcj) 



^0 



i / dtm^ie'"''!*! i / dt2n^'-'e"'-''^*^ ..A 
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(ki) . . . A^^ (km) • 



(5.19) 



Carrying out the U integrals to the right of the NE vertex, the coefficient of the product 
of gluon fields for this term reads 

Yl — — / diie^" *^i*i . . . / dtjtj -j- e'^-f^^*^ . 

(5.20) 

The right-most integral in this expression may be performed through integration-by-parts, 
and we find 



^NE,(3j) V'''' '^t) 
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This procedure can be iterated for the remaining integrations, and we find the expression 



ju2 ^fJ-i 

E__i_^_TT 
7(t, . KA^ 11 



2(n • K.)^ ^^n-Kk ^ 2(n ■ Kif ^^n-Kk 



(5.22) 



The first term in eq. (5.22) has the form of a product of eikonal emissions times a NE 
Feynman rule for gluon j, as expected. The second term implements correlations between 
gluon j and all gluons to the left of it. We may conveniently rewrite this to also involve all 
gluons to the right of gluon j, using the fact that one must ultimately sum over all gluons 
along the line. One may then use the identity 



Z^Z^ 0(„ . KA'i 11 
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(5.23) 



-43 - 



where we relabelled j -H- i in the second line. In this way we collect terms from all diagrams 
where j is either a single NE emission, or is correlated with gluons to its right. We can 
then replace eq. (5.22) with 



n 



n 



ri' 



W^? n-K, ^ 2(n • Kj)^ J-^. n ■ K, 



(5.24) 



This must be combined with the contributions arising from the second term in eq. (5.18), 
which generates diagrams of the form shown in fig. 13. This term correlates gluons at 
different positions along the line and, as before, the paired gluons are not necessarily 
adjacent, but may have additional eikonal emissions in between. Considering diagrams 
where the left-most gluon of the pair is the j^^, and performing manipulations similar to 
those leading to eq. (5.24), we are led to 



X - ki ■ kj n^J -i-r n 



\ ^ \ ^ rvi 



i>j 



ki ■ kj n'^' -i-r 

' n 



(5.25) 



The structure of this result is analagous to eq. (5.22), and follows from iterated integrations 
by parts. The second term implements correlations involving one gluon to the left of gluon 
j (gluon I), and one to the right of it (gluon i). 

As before, it is convenient to regroup terms in order to consider all diagrams where 
gluon j is the left-most one. One may do this after summation over all gluons j using the 
identity 

ki ■ kj n^^' \ - >. - ki ■ ki n^^ 



^ ^ in. Kjf 



(5.26) 



3 Kj i>j j j<l<i 

where we have interchanged the order of summation over j and I, and then relabelled j I. 
We can now replace eq. (5.25) with 



"NE, 
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(5.27) 



The total contribution from the terms of eq. (5.18), for each gluon j and its correlations 
with those to its right, as obtained from the sum of eqs. (5.24) and (5.25), is given by 



^ \ i>j j<l<i I k^j 
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iJ-k 



n ■ Kk 
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(5.28) 



This term has the form of a product of eikonal emissions, multiplied times a NE emission 
vertex for the emission of gluon j. This vertex has precisely the form of the second term 
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Figure 14: A diagram composed of two eikonal webs Wi and W2, correlated by a next-to-eikonal 
two-gluon vertex. Such diagrams contribute to the NE exponent, and thus can be described as NE 
webs. 

in the effective NE Feynman rule in eq. (4.17). This is a non-trivial check on the structure 
of this term, since the various scalar products ki ■ ki coming from eq. (5.18) must combine 
to produce precisely the squared partial momentum sum Kj. 

The last term to be considered is the second term on the second line in the exponent 
of eq. (5.6). In the non-abelian theory, and for A = 1, it gives 



The delta function implies that the two gluons are emitted from the same position on 
the eikonal line. Thus, this term correlates adjacent gluons, and indeed is easily seen to 
generate the NE two-gluon vertex Feynman rule of eq. (4.19). 

For gluons in the same group, we have now shown that the NE Feynman rules derived 
in this paper precisely correspond to those that arise from the path integral formalism of 
Ref. [39]. To complete the demonstration of this correspondence, some remarks arc in order 
concerning what happens when gluons belong to different groups. We will sec that this 
analysis gives a useful interpretation of the remainder term found in eqs. (3.49) and (3.55). 

In the path integral formalism, it can be proven that diagrams that are two-eikonal- 
line reducible do not contribute to the exponent. This is well-known for eikonal webs 
[11, 12, 10]: the non-trivial statement proven here and in [39] is that the same is also true 
for next-to-eikonal webs. One must however keep in mind that there are NE webs arising 
from diagrams such as the one shown in fig. 14, consisting of what would be a product 
of two eikonal webs at eikonal order, but where two gluons (one from each web) become 
correlated by a NE two-gluon vertex acting at different points on the external line. These 
diagrams are then two-eikonal-line irreducible, and thus must contribute to the exponent. 
Corresponding diagrams where the two-gluon vertex correlates gluons within one of the 
two webs are not present in the exponent, since they remain two-eikonal-line reducible. As 
a consequence, the contributions from all possible two-gluon vertex insertions no longer 
combine to form vertices depending on a single partial momentum sum, as they did above 




Jo Jo 




(5.29) 
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in the case of diagrams containing a single web. Instead, the leftover contributions form a 
remainder term, which has the same form as that discussed previously, in Sec. 3.4. 

In more detail, consider again eq. (5.17), and assume for simplicity that the line being 
considered has contributions from two groups g and h only (the following discussion easily 
generalises to the case of more groups). Let us also take gluon j as coming from group g. 
One may then split the sum on the left-hand side of eq. (5.17) into two parts, involving 
gluons from g and h respectively, to get 

i>j,i&g i>j,i&h I k^j ^ ' 

where as usual gi is the momentum of the i*'' gluon from g. The first two terms in brackets 
combine to give the combination of partial momentum sums Gj+Gj+i, as already discussed 
in eq. (5.30), and lead to a next-to-eikonal emission vertex for gluon j with eikonal Feynman 
rules for the remaining gluons. The second term can be rewritten as 

-(W+G,)liA"^J' ^'-''^ 

where i labels the gluon from h that is nearest (on the right-hand side) to the gluon 
from g. This has the form of a two-gluon correlation, with eikonal Feynman rules for the 
remaining gluons, and is the same as that given in eq. (3.49). Note that the previous result 
is expressed as a sum over all possible permutations of the gluons along the line, whereas 
here we consider only one such permutation. Also, we have yet to sum over the gluons j 
from group g. 

One may also consider eq. (5.28), again in the case in which two groups are present 
along the line. Considering the bracketed combination in the first line, all terms involving 
gluons from the same group combine to form a single squared partial momentum sum 
or Hj, as already discussed. There remain the terms in the double sum, in which the two 
gluons come from different groups. Considering the case where gluon j comes from group 
g, the leftover terms are 

m 

2Y,Y.9rhk = 2Gj-Hi, (5.32) 

i=j k>i 

where again i labels the gluon from h that is nearest (on the right-hand side) to the gluon 
from g. Combining with the other factors in eq. (5.28), one has a contribution to the 
remainder term of 

Gj ■ 11, TT "''M 
2n • (G,- + H,) ii ■ ^^-^^^ 

After summing over all gluons from 5, and all permutations of gluons on the line, this 
agrees with the result obtained in eq. (3.55). We have therefore demonstrated the full 
equivalence between the results of this paper and the path integral approach of Ref. [39]. 

The above discussion can be summarised as follows. In the diagrammatic approach of 
this paper, we find that contributions from different groups (derived from webs at eikonal 
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Figure 15: Born level diagram for the Drell-Yan production of a an off-shell photon. 



and NE order containing local vertices only) factorise, up to a remainder term that corre- 
lates gluons between different groups. The NE Feynman rules thus obtained are precisely 
those that one finds through the path integral formalism, where the remainder term can 
be interpreted as correlating pairs of eikonal webs to form a new class of NE web. That 
the same structure and effective Feynman rules emerge from both the diagrammatic and 
path-integral viewpoints, is a highly non-trivial check of both approaches. 

6. An illustrative example: the Drell-Yan process 

In order to illustrate the results obtained in the previous sections, showing how they can 
be implemented in a concrete calculation, we consider in this section soft gluon corrections 
to the Drcll-Yan production of a virtual photon in quark-antiquark annihilation. Since 
our focus here is just to give an example of the implementation of our NE rules, we will 
restrict our attention to abelian-like contributions to the cross section, which build up 
the color structure at N^LO. Furthermore, we will consider only real radiation, since 
when including mixed real-virtual corrections we will need to account for Low's theorem 
[41, 42, 43] contributions, which are not within the scope of the present paper. With these 
restrictions, wc will compute soft gluon corrections up to NE level, and up to NNLO, 
using the rules developed in the previous sections. This is the simplest situation in which 
one can explicitly test NE exponentiation (as some of the two-loop NE corrections arise by 
expanding the exponential of the NE one-loop result) . Furthermore, this calculation can be 
viewed as a preliminary step towards the generalisation of existing soft-gluon resummation 
methods to NLO in the threshold expansion, as well as a check of the validity of our 
analysis. 

The NLO corrections to the Drell-Yan process were computed in [48], while the NNLO 
results was obtained in [49] (eikonal NNLO results were first presented in [50, 51]). We 
begin with a very brief summary of the LO calculation in order to introduce notation. 

The LO Drell-Yan process [52], shown at parton level in fig. 15, proceeds through the 
channel 

Q(p) + Q(p)^7*('?"), (6.1) 

where Q denotes a quark, and is the momentum of the off-shell photon. Since at higher 
orders (or at hadron level) the ccntcr-of-mass energy s and the photon mass = do 
not coincide, it is customary to parametrize the Drell-Yan total cross section in terms of 
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the variable z = Q'^/s,Q<z<l. One then normaUscs the higher order cross-section to 
the Born result defining the n-loop K factor X*-"^ (z) as 

(6 2) 

where C7(") is the cross-section at 0{ag). The Born amplitude is simply 

■^fo) = u{p)rv{p) , (6.3) 

where we drop coupling factors and spin averaging, as these will ultimately cancel upon 
taking ratios to compute the K factor. Up to inessential factors, the Born cross section is 
given by 

jd. 



/ (2^ ^^^^ {<1-P-P)S+ {q' - Q') TV [|)7"^7a] 

= 2'K5+{q^ -Q^)Ap-p{d-2), (6.4) 

where we will henceforth take d = 4 — 2e. Clearly, the K factor at leading order is simply 
given by 

K^^\z) = 5{\ - z) . (6.5) 
We now move on to the NLO, next-to eikonal calculation. 

6.1 Next-to-eikonal contributions at NLO 

So far in this paper, we have only considered NE contributions to matrix elements. In 
calculating soft gluon contributions to cross-sections, one must also consider phase space 
integration. In the eikonal approximation (as is well known), the n-gluon phase space can 
be written as a product of decoupled one-gluon phase spaces. This is a necessary condition 
for the exponentiation of the cross-section which underlies soft-gluon resummations. At 
NE order the exact factorization of phase space breaks down, and it becomes necessary 
to study the unfactorized corrections to all orders. This is done, for the inclusive cross 
section, in Appendix B. Preliminarily, we note that at NE level one expects the schematic 
structure 

(TNE = / dPSeik IMneI' + / dPSNElMeiki' • (6.6) 

Here A4eik and A^ne are the eikonal and next-to-eikonal contributions to the matrix ele- 
ment, extensively discussed in the previous sections, and we have introduced an analogous 
expansion for the multi-gluon phase space measure of integration dPS. The NE correction 
to the phase space measure will also induce correlations between pairs of soft gluons, with 
a sum over all possible correlations. Up to NNLO, these corrections are easy to implement 
since the exact phase-space integration can be performed analytically, as we will do below. 

We begin by computing the eikonal contributions to the NLO ii'-factor. Eikonal Feyn- 
man rules easily yield 

= 2glCF—{\-e). (6.7) 



•^eik,a 
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Figure 16: Squared matrix elements contributing at NLO in the eikonal limit. The cut denotes 
the final state. 



where the factor of 2 in the first line comes from combining the two diagrams shown in 
fig. 16, and we have introduced the Mandelstam invariants 

t = -2k-p; u = -2k-p. (6.8) 

To carry out the phase space integration, one may parametrise the four-vectors as 



P 



^(1,1,0^); p=^(l,-l,0^) ; A: = -L^^(l,cos^,sin^,0). (6.9) 



Defining then 



one finds 



1 + cos ( 



(6.10) 



t = -2s{l-y){l- z); u = -2sy{l - z) . (6.11) 

The NLO i^-factor, defined in cq. (6.2), is computed by including a factor 5{k^ ~ (1 ~ z)p^) 
in the two-particle phase space, and introducing the appropriate normalization, according 
to eq. (6.4). Setting the renormalization scale /x^ = one finds, at the eikonal level 



M 



(l)a 
eik 



(6.12) 



Substituting the squared eikonal matrix element given in eq. (6.7), performing the integra- 
tion and Taylor expanding in e this yields 



a 



4k (^) = ^Cf{- -Vo{z) + 16Pi(z) 



47r 

•4e 



81og(2:) 
1-z 



AT, ( \ -i/- ( \ 4l0gzl0g(l - Z) log2 Z 



(1-z) 



1-z 



(6.13) 



where, following [49] , we use the notation 

Vp{z) = 



" log^(l-^) " 
l-z 



(6.14) 
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Figure 17: Diagrams contributing at NE level to the squared amplitude for NLO real emission. 
Two additional diagrams, arising from complex conjugation, are not shown. 



to denote plus distributions. Note that we are not displaying terms proportional to (5(1 — z), 
which mix with virtual contributions not considered here. Note also that eq. (6.13) contains 



terms that do not behave as (1 



-1. 



since we are using eikonal matrix elements, such 



terms arise from subleading corrections to the phase space measure. 

Let's now consider NE corrections to the matrix element. According to the results of 
Sec. 3, the effective Feynman rule for single gluon emission at NE order is 



2p-k 



(6.15) 



For real emission, the second term in eq. (6.15) vanishes since k'^ = 0. The diagrams 
contributing at NE level arc shown in fig. 17, and consist of all possible insertions of at 
most one NE vertex in the eikonal diagrams of figure 16. The corresponding contribution 
to the squared matrix element is 



M 



(1) 

NE,a 



M 



eik 



t 



(1) 



eik,a 



M 



(l)a 
NE 



t 



p ■ k 
8(1- 



2p-k 



p ■ k 2p ■ k 



(6.16) 



Carrying out the phase space integration, dividing by the Born normalisation and expand- 
ing in e leads to 



a. 



— 16 log(l — z) + S log z 



-4e [-41og^(l - z) +41og2;log(l - z) - log^ z + 3(2] 



(6.17) 



Combining eqs. (6.13) and (6.17) and retaining only terms with a logarithmic dependence 
on z, up to NE order, we obtain 

= ^Cf \ - -Vq{z) + 16Pi(z) - 161og(l -z)- ieUv2iz) - 3C2Po(^) 
47r I e L 



-41og2(l-z) + 41og(l-z) 



(6.18) 



Note that to obtain this result we expanded around z = 1, including for example the 
replacement log z — >■ — (1 — z). 

One may compare this result with the corresponding exact NLO result including the 
complete z dependence [48]. One easily verifies that expanding the exact result around 
z = 1 gives eq. (6.18). Thus, the effective Feynman rules derived in the previous sections 
pass a first simple test. 



6.2 Next-to-eikonal contributions at NNLO 

In this section we compute the contribution proportional to Cp to the NNLO if-factor 
for DY production of a virtual photon, and, as in the previous section, we restrict our 
attention to real gluon emission from a qq initial state. This is the simplest possible 
example illustrating the analysis performed in sections 2-4. 

The structure of NE corrections found in the preceding sections implies that the am- 
plitude for the double-real-emission graphs can be extracted, up to NE accuracy, from the 
expansion of an exponential with the schematic structure 



7W = A^(o) exp 



,(2) 



(6.19) 



where M.^^^ oc Tr [^7" pja]^ £^nd the bar over an amplitude represents the expression for 
that amplitude with Born result factored out. Note there is no O(a^) eikonal amplitude 
in the exponent. This is because we are considering only abelian graphs, and in that set 
there is no connected subdiagram involving two eikonal real emissions. In this language, 
the single and double real radiation matrix elements, up to NE accuracy, can be written as 



M 



(1) 

eik 



M 



NE 



NE 



(6.20) 



Since wc arc not considering virtual corrections, the only contributions to the squared 
matrix element that we will need at 0{a1) and at NE accuracy are given by 



At (2) J^{2) 



t 



-F2Re 



-F2Re 



M 



(2) 
NE 



M 



ik-'^'NE J 



eik 



+ 



NE,lg 



+ 



NE,2g 



(6.21) 



where, as the notation suggests, the first term is the eikonal approximation to the double 
real radiation matrix element squared, the second term is the NE correction arising from 
single-gluon NE vertices, and the third term is the NE correction arising from double-gluon 
NE vertices. We will treat each contribution in turn. 

Let's begin with the eikonal result, which has the explicit expression 

M'^^^^ = Af — ; ^''^ , _ . TV [^7" ^7«] , (6.22) 



eik 



p ■ kip ■ k2P • kip ■ k2 
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where, for simplicity, we have included in the normalization N all coupling, color and spin 
averaging factors, most of which drop out when constructing the K factor. 

The integration of eq. (6.22) over the three-body phase space is non-trivial, and we 
briefly discuss it here, since it will be needed also for the evaluaton of the NE corrections. 
We follow Ref. [49], and we evaluate the integrals in the center-of-mass frame of the two 
outgoing gluons. In order to perform the integral, we parametrise the momenta in d = 4— 2e 
dimensions as in [49]. We write 



k2 
P 

Q 

P 



2 

^/sT2 
2 

is - i) 

2-^512 



(1,0, . . . ,sin^2sin^i,cos02sin6'i,cos^i) , 
(1, 0, . . . , — sin02 sin^i, — cos 02 sin^i, — cos ^i) 
(1,0,... ,0,1) , 



(6.23) 



Sl2 



2^/^ 



12 



,0, ... ,0, jgjsinV', IgjcosV' , 



i + si2 - 



2^/iIi 

Here we have introduced the quantities 



,0,...,0, IglsinV', l^cos^;- 



2^/iIi 



i=2p-q = ip + qf-Q\ 
u = 2p ■ q = {p + qf - Q"^ , 
si2 = 2ki ■ k2 = s — i—u + Q'^, 
{s - Q^){u - Q^) - suit + 



(6.24) 



COSl/j = 

k1 = 



(s-t)VA(s,Q2,si2) 
\/X{s,Q^,si2) 



where A is the Kallen function, A(a, b, c) = + 6"^ + (P — 2ab — 2ac — 26c. The Mandelstam 
invariants t and u can in turn be expressed as functions of the photon energy fraction 
z = Q^/s and of two further variables x and y, normalized so that < a; < 1 and 
< y < 1, with the definitions 



s[l-y{l-z)] , 



t = s 



z + y{l- z)- 



y{l - y)x{l - z) 



l-y{l-z) 

The d-dimensional three-body phase space can then be written as 



(6.25) 



dPS. 



-il-x) 



2d-b 



Jo Jo 



(102 (smOi)'^-^ (sin 02 



(47r)^ r(d - 3) 

X [\y I'dx [y{l - y)] [x{l - x)] [l - y(l - z)\ 
Jo Jo 



\d-4 



l-d/2 



(6.26) 
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To integrate the matrix element squared in eq. (6.22) we still need expressions for the 
propagator factors in the current parametrization. They are 



P ■ ki 

P ■ k2 
p-ki 
p-k2 



s-t 
4 

s-i 



(1 - cos Oi) 
(1 + COS 6*1) 



A — B cos 61 — C sin 9i cos 62 
A + B cos 9i + C sin 9i cos 62 , 



(6.27) 



where 



A 
B 
C 



t + s 



12 



I q\ cos ^p — 



/Sl2 



\q\ smi/j . 



(6.28) 



It can be verified that = B^ + C^. Using this fact and putting together eq. (6.22) and 
eq. (6.26) one may verify that all the angular integrals to be performed can be reduced to 
the standard form 



sin*^ ^ 9i sm' 



d-^e2 



cos di)P{l — COS X COS 9i — sin x sin 9i cos ^2)^ 



(6.29) 



r(d-2-p-g)r2(f-i) 



d , 2(X\ 



where 2-P1 is the standard hyper geometric function. Identifying cosx = —B/A, using 
cos^(x/2) = (1 + cosx)/2, and expanding in powers of e = 2 — d/2, one gets, for example, 



7-(l,l) 

'4-2€ 



(x) 



TT 

e 



A + B 
2A 



-l-e 



1 + Li2 



A-B 
2A 



(6.30) 



After performing the angular integrals in this way, the remaining integrations over x 
and y are easily carried out, after expanding the integrand in powers of (1 — z). This is 
consistent here since, for real emission contributions to the Drell-Yan cross section, the soft 
gluon expansion and the threshold expansion coincide. The resulting eikonal contribution 
to the NNLO K factor reads 



32 ^ , ^ 128 ^ , , 256 ^ , , 320 , 
-^Vo{z) + —Vi{z) -—V2{z) - — log(l 



+ ^ V^{z) + 6401og2(l -z) + ... 



(6.31) 



Note that for the sake of illustration we display only terms with a logarithmic dependence 
on (1 — z), and specifically those with rational coefficients, that are easily extracted from 
the exact calculation for comparison. 
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fa) (b) 



Figure 18: Cut Feynman diagrams involving the next-to-eikonal two-gluon vertices summarized 
by the tensor R^"^. Note that also the complex conjugate diagrams must be included. 



We now turn to NE corrections to the matrix element, and specifically we consider the 
second term in the last line of eq. (6.21), which arises from the insertion of single-gluon 
NE vertices. The application of the rules derived in the previous sections yields 



NE,lg 



+ 



+ {ki -H- k2) 



1^[^7"^7a]. (6.32) 



_p- kip- ki \p- k2 p- k2. 

The integration over phase space can be carried out as for the eikonal term, and one finds 
a contribution to the K factor given by 



(2) 

NE,lg 



1024 



12 



256 



log(l -z) log {1- z 



+ — — log'^(l-2;) + ... 



(6.33) 



where again for simplicity we display only logarithmic terms with rational coefficients. 

Finally, we need the contribution given by the final term in eq. (6.21), which involves 
NE contributions to two-gluon emission that don't already arise from the exponentiation 
of the single emission. That is, one must include the contribution of the two-gluon vertices 
discussed in Sec. 3. The sum of the relevant vertices on, say, the quark line carrying 
momentum p gives the tensor 

{p ■ k2)p^kl + [p ■ ki)k^p'' - {p ■ ki){p ■ k2) g^'' - {ki ■ k2)p^p'' 



R^'^kiM-.P) 



(6.34) 



p- {ki + k2) 

where an additional factor of two arises upon summing over the different orderings of the 
gluon momenta. Note that, given that this vertex joins two adjacent gluons, it includes 
the seagull vertex (i.e. the term in g^^), as discussed in section 4.1. It is easily checked 
that this gives zero if either gluon lands on a fermion line carrying the same momentum 
p as the emitter. Thus, the only possible contributions come from diagrams such as those 
shown in fig. 18. In fact, these diagrams also give a vanishing contribution, as can be seen 
from the following argument. Consider for example the diagram in fig. 18(a). This involves 
the factor 

R^^PtiPu P ■ P 



p ■ kip ■ k2P ■ kip ■ k2 



p- {ki + k2) p ■ kip ■ k2P ■ kip ■ k2 

X p ■ k2P ■ ki + p ■ k2P • ki — ki ■ k2P ■ p 



.35) 
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Note that the term in i?^^ involving g^^ vanishes, since = 0. If we now pick hght-cone 
coordinates, such that = {p^ ,0~ ,0±) and p^ = [0'^ ,p~ ,0±), with p^ = [p^ ±p^)/\/2, 
we observe that the factor in brackets on the r.h.s. of eq. (6.35) becomes 

g 

p ■ k2P ■ h + p ■ k2P ■ h - ki ■ k2P ■ p = - 2 '^i.-L ■ ^2,± • (6.36) 

The prefactor of this term depends only on longitudinal momentum components. As a 
consequence, the diagram involving the two-gluon vertex gives a vanishing contribution 
upon integration over the gluon transverse momenta. 

We conclude that the NE contribution to the NNLO K factor is given just by the sum 
of eqs. (6.31) and (6.33). This, as expected, agrees with the result obtained from a full 
calculation of the double real emission contribution, followed by an expansion to NE order. 
The latter result is not available directly in the literature, but we include a short summary 
of the calculation in appendix C. 

7. Conclusion 

In this paper we have considered the exponentiation of soft gluon contributions to gauge 
theory amplitudes and cross sections, from a Fcynman diagram viewpoint, classifying the 
structure of ncxt-to-cikonal corrections. Wc showed that the notion of exponentiation can 
indeed be extended to NE order, as given by eq. (4.16). We note however that at NE level 
there are non-factorizable contributions which must be separately considered, and which 
are related to the radiationless amplitude by Low's theorem and its generalizations. Fac- 
torizable contributions exponentiate and fall into two classes. First, there are conventional 
(eikonal) webs, modified by including effective NE Feyman rules; then, there are new webs 
formed by taking pairs of eikonal webs, and correlating a pair of gluons (one from each 
original web) with a two-gluon vertex. The derivation of this result relies upon establishing 
the appropriate generalisation of the eikonal identity, which is given in eq. (3.24) and ex- 
presses the fact that contributions from different groups factorise on a given external line, 
up to a remainder term which gives rise to the two-gluon vertex contributions mentioned 
above. 

We performed a detailed comparison of the effective NE Feynman rules derived in this 
paper with the results obtained using the path integral formalism of [39] , finding complete 
agreement between the two formalisms. This is a highly non-trivial cross-check of our 
classification of NE contributions to matrix elements, and also provides further insight on 
how to apply them in practice. In particular, a diagrammatic formalism makes contact 
with conventional factorization proofs, and is likely to be amenable to generalizations to 
hard coUinear emissions which would be needed in order to perform a complete analysis of 
the threshold expansion beyond leading order. 

The complete structure of next-to-eikonal corrections to QCD cross sections receives 
contributions not only from matrix elements but also from the integration over phase space 
for multiple soft-gluon emission. The discussion in this paper (see appendix B) shows how 
to treat these corrections for any number of gluons in the case of the Drell-Yan process; 
we expect that it should be possible to generalize these arguments to other processes. 
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Having classified the structure of NE corrections, it is interesting to ponder what might 
happen at NNE order (and beyond). Terms at this order in the expansion in gluon momenta 
are no longer divergent, but may nevertheless be numerically significant. The path-integral 
arguments of [39], as well as the recursive arguments presented in this paper, make clear 
that there is always a subset of external emission graphs which exponentiates at any given 
order in the momentum expansion. At higher orders, new vertices involving multiple gluon 
emissions will arise. More serious problems are bound to emerge, however, from non- 
factorizable contributions, which are governed by Low's theorem at NE level. It is difficult 
to generalise the Low-Burnett-KroU argument to higher orders in momentum, as gauge 
invariance does not then completely fix the structure of non-exponentiating contributions. 
Nevertheless, it might be the case that factorizable contributions dominate in practice, 
allowing for at least a partial resummation of NNE terms and beyond. It should also be 
noted that, as more finite orders terms are generated through exponentiation, it would 
become increasingly important to accurately match the resummed result to fixed-order 
perturbation theory, since a simple expansion of the exponential is bound to generate 
an increasing number of subleading terms which are beyond the accuracy guaranteed by 
resummation. 

The results of this paper should be seen as a step towards the generalization of Sudakov 
resummation beyond leading order in the threshold expansion. In order to achieve the goal 
several further steps are still needed: the simple organization of phase space proposed 
in this paper must be brought to bear in a more formal way and in more general cross 
sections; non-factorizable contributions must be better understood; finally modifications 
of collinear evolution for hard collinear emissions must be fully understood at next-to- 
leading order in the threshold expansion. In order to show that our work is a step in 
the right direction, and as a first illustration of how to apply our effective Feynman rules 
at NE order, we reproduced a subset of NE logarithms in the real emission contribution 
to the Drell-Yan K-factor up to NNLO (specifically, those associated with the Cp color 
structures). This is in a sense the simplest possible example, given that the limit z — >■ 1 
for this process corresponds to soft singularities only, with no mixing from hard collinear 
divergences. One must also note that we have only discussed the simple case of a color- 
singlet hard interaction linking two partonic lines. We expect however that the bulk of our 
considerations will extend to more general color configurations, as indicated in [13]. Work 
towards understanding these various generalizations is ongoing. 
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A. Rearranging NE numerators for spin-^ particles 

We wish to prove that, when acting on a massless spinor, the following identity holds: 



jyfj,l...fj.n ^ 

1=1 

n 

= E 



1=1 



(A.l) 



where we use N^^^-^^" as a shorthand for —J\f^^"2'^{jp,ki), and where Ki = Ylp=i^p- ^otc 
that, as the notation suggests, the second term in the second line vanishes for i = 1. We 
begin by observing that cq. (A.l) can easily be verified for low values of n, and indeed is 
trivial for n = 1. This suggests setting up a recursive argument. In order to proceed, we 
first isolate the contribution of the first gluon in the l.h.s. of eq. (A.l), writing 



J=2 



i=2 



i-l 



j=i+i 



k=2 



(A.2) 



We can now use the induction hypothesis, assuming that eq. (A.l) holds forn — 1 emissions, 
and using it to replace the sum in the second term of eq. (A.2). We get 



(A.3) 



j=2 



i=2 



Observe that the factor {2pi^^) never appears on the second line, where the sum starts 
with the second emission, and for the same reason the last term starts contributing only 
for i = 3. We now operate on the first term on the second line of eq. (A.3), using the 
commutation relation 

[pl^, = 2p-k -ff^Y -2p'' 7^^^ (A.4) 

in order to bring the operator to the right of the expression, where it annihilates the 
massless spinor. We find 



j=2 



1=3 



n i2pl'i)2p■Ki-^''^-^^^'^ 



(A.5) 
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The last term in the sum on the first hne has a factor 0^=2 (2p^^ ) for every value of the 
summation index i: one can then perform the sum getting simply a factor of -^2- This 
can be combined with the first term in eq. (A. 5), using the fact that Ki — K2 = ki. The 
result is simply 11^=2 {'^P^^)l^^^i- This can, in turn, be combined with the first term in 
the sum on the first line, where the factor of p can be commuted across 7'^^ to yield 2p'^i . 
One now has precisely all terms needed to reconstruct the sum of the first terms on the 
r.h.s. of eq. (A.l). One finds at this stage 



E 

1=1 



\{{2p 



i=2 



J|(2p'^^) 2p- fen"'?"* 



i=3 



(A.6) 



Next, we act on the second line in eq. (A.6) applying the commutation relation 

[^'', = Yi^ '^f - Yi"" 2 p^" , (A.7) 
also in order to bring the factor to act on the massless spinor to the right. This leads to 



i=l 

n 

+ E 



i=2 



J](2p^0 2p-ka'''7''' 



(A.8) 



1=3 



At this stage, one just needs to properly recombine the terms in in the second and in the 
third sum in eq. (A.8) in order to recover eq. (A.l). We do this in the following way. 

• Combining the term with i = 2 in the second sum with the term i = 3 for the second 
term of the third sum, and using the fact that k2 + K3 = K2, one gets from these 
two terms a contribution 



J](2p'^^) 2p-K2Y'Y\ 



(A.9) 



which is the first non-vanishing contribution (the one with i = 2) to the last term in 
eq. (A.l). 

• One easily recognizes that all remaining terms in the last sum in eq. (A.l) are given 
precisely by the first term in the last sum in eq. (A.8), when the sum is performed 
in the range 3 < i < n), as written. Eq. (A.l) is thus already complete: it remains 
to be shown that the leftover terms in eq. (A.8) cancel. 

• Consider indeed the remaining terms in the second sum in eq. (A.8) (those with 
3 < i < n). Using the fact that ki — Ki = — Kj+i, they can be combined with the 
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third term of the third sum, which is summed over the same range, to get 



i=3 



(A.IO) 



It is easy to see that this sum precisely cancels the leftover contributions from second 
term in the last sum in eq. (A.8), which are in the range 4 < i < n, by just shifting 
i + 1 ^ i in eq. (A.IO). This completes the proof of eq. (A.l). 



B. Next-to-eikonal phase space at arbitrary order 

In the Sec. 6.2 we have illustrated the use of next-to-eikonal Feynman rules by calculating 
certain NNLO contributions to the Drell-Yan cross section. In order to correctly calculate 
logarithmic contributions which become dominant near partonic threshold, one must in- 
clude corrections to the phase space measure of integration, as well as to matrix elements. 
We have seen that, at one and two loops, these corrections have a simple form and are eas- 
ily taken into account. In order to perform a resummation, however, we will need to study 
NE corrections to the phase space measure to all orders in perturbation theory, that is for 
any number of emitted gluons. In this appendix we show, for the Drell-Yan process and 
for obscrvablcs where the soft gluon phase space is integrated over, that the phase space 
measure remains very simple at NE order and for any number of emitted gluons. Indeed, it 
turns out to be possible to express this measure assigning identical factors to each emitted 
gluon, in a way that generalizes the known factorization properties of multi-gluon phase 
space at eikonal level. This will be instrumental in extending resummation techniques for 
physical cross sections to subleading orders in the threshold expansion. 

Let us briefly recall the effects of subleading phase space corrections at one and two 
loops. In the case of one final state gluon (eq. (6.12)), dependence on the energy fraction 
z factors out of the integration over the angular parameter y, and subleading corrections 
straightforwardly result upon expanding the z-dependent prefactor in powers of e. When 
two real gluon emissions are present, as in eq. (6.26), the integration measure contains the 
factor 

[l-y(l-z)]-'+^ = l + (l-e)(l-z)+0[(l-z)2] . (B.l) 

Keeping only the first order term on the right-hand side corresponds to the NE approxi- 
mation, which in the present case coincides with the next-to-lcading order in the threshold 
expansion. We will now show how eq. (B.l) can be simply generalized to n-gluon emission. 

Our starting point is the (n-|- l)-body phase space for the emission of n gluons together 
with a final state vector boson. Using standard factorization properties, we write this as 

f 1'°^ dK"^ I' 

J dPS{p + p^q + ki + ...kn) = J J dPS{p + p-^q + K) 

X J dPS{K ^ ki + ...kn) . (B.2) 

Here dPS(^ — > B) denotes the phase space for the reaction A ^ B, q is the 4-momentum 
of the vector boson and {ki} the momenta of the emitted gluons. On the right-hand side of 
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eq. (B.2) wc have introduced the vector sum of tlie emitted gluon momenta K = X^iLi 
One may interpret the product of phase spaces on the right-hand side of eq. (B.2) in terms 
of the production of a pseudo-particle of invariant mass JC^, followed by its decay into n 
massless gluons; one then integrates over all possible invariant masses K"^. The convenience 
of this decomposition is that the two body phase space in eq. (B.2) is the same for any 
number of gluons, and depends upon the momenta of the latter solely through K^. One 
may then evaluate the two-body phase space by parametrising momenta as 

p= ^(1,0,... ,0,1) ; 

^= ^(1,0,..., 0,-1) ; (B.3) 



? = (vl9l^ + Q^,0, •••J^sinV', IglcosV'^ 
K= (VkT + 0, . . . , -Igl sin V', -k1 cos V') . 
Straightforward manipulations in d = 4 — 2e dimensions yield 

where A(a;, y, z) = + + — 2xz — 2xy — 2yz is the Kallen function. As in the case 
of one real gluon emission (corresponding to = 0), one may define the Mandelstam 
invariants 



t = -2K-P = -^/i(VkT + Q'-|^cosVj ; 

u = -2K-P = -V^(Vl^^ + <3^ + l^cosV) . (B.5) 
One may then show that 



2 , 4{tu-Q'^s) 



so that the two-body phase space may be written as 

dPS{p + p^q + K) = s-^+^ [tu - Qh) '/'-^ di, . (B.7) 

The n-body phase space in eq. (B.2), on the other hand, is complicated in general. If 
however one performs all integrations for fixed total invariant mass K^, the dependence of 
the result on is dictated by dimensional analysis, and is of the form 

j dPS{K^ki + ... kn) = (i^2^"-2-(n-l)e j dPS {K ^ ki + . . . kn) , (B.8) 

where dPS on the right-hand side denotes a dimensionless bounded measure of integration 
depending only on the angles between the gluons in a chosen frame (for example their 
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center of mass frame, as in Sec. 6.2). This measure is therefore independent of K^, and 
z = Q'^ I s. Using the Jacobian 

dK^di\) = ^ dtdu, (B.9) 
ytu — Q^s 

the full phase space of eq. (B.2) becomes 

JdPSip + P^q + h + ...kn) = j^ S-'+' J dt du (i^2)n-2-(n-l)e 

X {tu - Q^sY' J dPS{K ^ ki + ...kn) . (B.IO) 

As in Sec. 6.2, it is convenient to introduce the variables x and y defined implicitly in terms 
of u and t as in eq. (6.25). The Jacobian is 

dudt = yC^-y) (i-zfs'^dxdy, (B.ll) 
l-y{l-z) 

so that the full phase space takes the final form 



/ 



dPS(p + p^g + fci + ...M = ^1^2 r|i?e) ^""'""(1 - zf^-^-^^' (B.12) 



/' dy C dx [y{l - y)^-'-'^' x'n-^-in-De (1 _ ^ye 
Jo Jo 

[1 - y{l - z)] JdPS{K^ki + ...kn) . 



One sees that the dependence on the energy fraction z has only two sources: the over- 
all prefactor of (1 — 2;)^"'~^~^"'^, which conspires with the matrix element to give the 
leading power (1 — z)~^ in the threshold expansion of the cross section, and the factor 
[1 — y(l — z)]~^"~^''^^~^\ which is neglected in the eikonal approximation. One observes 
that eq. (B.12), for n = 2, reduces as expected to the three-particle phase space measure 
in eq. (6.26). 

It is interesting to note that, at NLO in the threshold expansion, the z dependence of 
the phase space measure in eq. (B.12) can be expressed in factorized form, with the same 
factor assigned to each emitted gluon. Indeed 

[1 - y{l - z)] = 1 + _ 1) y (1 _ ^) (1 _ e) + o[{l - zf] 

= z^^-'^y (^z-^^-'^^y . (B.13) 

This provides, to some extent, a generalization of the factorization property of the multi- 
gluon phase space, which is known to degenerate into a product of uncorrelated single- 
gluon phase spaces at leading order in the threshold expansion. At NLO in the threshold 
expansion, the multi-gluon phase space measure is still factorizable, however the factors 
assigned to each emitted gluon retain a dependence on the angular variable y describing the 
overall orientation on the multi-gluon system. In other words, single-gluon phase spaces 
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Figure 19: Diagrams for tlic doublc-rcal-cmission contribution to the NNLO Droll- Yan K factor 
discussed in the text. A cut is implied over the intermediate state in each case, and complex 
conjugates of the above diagrams (excluding (e), which is real) must also be included. 



are not fully uncorrelated, but their correlation is simple, depending only on the global 
variables of the multi-gluon system and not on individual gluon momenta. 

This discussion applies to the explicit example of Drell-Yan production. We expect 
that such arguments will apply more generally in other scattering processes, pending a 
suitable parametrisation of the partonic momenta. 



C. The double-real-emission contribution to the Drell-Yan K factor 

In this appendix we briefly describe how to compute the terms proportional to C|, of the 
Drell-Yan if -factor, for the qq initial state, by using ordinary Feynman diagrams and ex- 
panding them to NE order. The relevant diagrams are shown in fig. 19. The corresponding 
squared matrix elements are easily computed, and must then be integrated with the phase 
space measure in eq. (6.26). As an example, diagram (a) contributes a factor 



1^,2 Tr[^7"(^- h- h)Y{^- h)l^ M- P+ h+ h)M- P+ h)lA 



Note that the contributions from diagrams (a) — [d) must be counted twice in order to 
include Hcrniitian conjugate graphs, while diagram (e) is real. 

To calculate the squared matrix element to NE order, one first relabels ki -> ^fcj, so 

that 

p-ki^^p-ki, p-ki^ ^p-ki, ki- kj ^'^ ki- kj . (C.2) 
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One then expands each diagram to first subleading order in ^, which corresponds to the 
NE approximation. Through repeated use of the identities 



p ■ ki 1 s — t 

I, (C.3) 



p- k2 p-k2 2 

p-ki 1 si2 + i-Q'^ 



p-k2 p-k2 2 

(with similar results for fei o A;2), each diagram can be written as a sum of terms containing 
no more than two factors of p ■ ki and p-ki. Then each term becomes an integral of the 
form of eq. (6.29). The remaining phase space integrals can be carried out after expanding 

the integrand in powers of 1 — z and e, as described for the NE calculation in Sec. 6.2. The 
final result for the full amplitude (keeping only logarithmic terms with rational coefficients 
as done in the text) is given by 



(2) . s / a.- ^ 2 



32 ^ , , 128 ^ , , 128 , 
■^Voiz) + —V^{z)--^ log(l - z) 



(C.4) 



256 , , 256 , 320, , 

- — V2{z) + — log2(l -Z)-— log(l - Z) 

+ ^ V^{z) - ^ log=^(l -z) + 6401og2(l - z) 
which is in complete agreement with the sum of eqs. (6.31) and (6.33). 
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